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ABSTRACT 

In solar euv spectra the He I and He n resonance lines show unusual behaviour and 
have anomalously high intensities compared with other transition region lines. The 
formation of the helium resonance lines is investigated through extensive non-LTE 
radiative transfer calculations. The model atmospheres of Vernazza, Avrett & Loeser 
(1981) are found to provide reasonable matches to the helium resonance line intensi- 
ties but significantly over-estimate the intensities of other transition region lines. New 
model atmospheres have been developed from emission measure distributions derived 
by Macpherson & Jordan (1999), which are consistent with SOHO observations of 
transition region lines other than those of helium. These models fail to reproduce the 
observed helium resonance line intensities by significant factors. The possibility that 
non-Maxwellian electron distributions in the transition region might lead to increased 
collisional excitation rates in the helium lines is studied. Collisional excitation and 
ionization rates are re-computed for distribution functions with power law suprather- 
mal tails which may form by the transport of fast electrons from high temperature 
regions. Enhancements of the helium resonance line intensities are found, but many of 
the predictions of the models regarding line ratios are inconsistent with observations. 
These results suggest that any such departures from Maxwellian electron distributions 
are not responsible for the helium resonance line intensities. 

Key words: line: formation - radiative transfer - Sun: transition region - Sun: UV 
radiation. 



1 INTRODUCTION 

The resonance lines of He I and He II show unusual behaviour 
when compared with other strong emission lines in solar 
euv spectra. Recent results from the Solar and Heliospheric 
Observatory (SOHO) show that their intensities in coronal 
holes are factors of 1.5-2.0 smaller than in the quiet Sun (Pe- 
ter 1999; Jordan, Macpherson & Smith 2001), while other 
lines formed at similar temperatures show only very small re- 
ductions in intensity. In the quiet Sun, the helium resonance 
line intensities are at least an order of magnitude too large 
to be reproduced by emission measure distributions that ac- 
count for other transition region (TR) lines (Macpherson & 
Jordan 1999 - hereafter MJ99). For discussions of earlier 
observations of the helium resonance lines and attempts to 
explain them see Hammer (1997) and MJ99. The discrepan- 
cies between observations and modelling of the helium lines 
in the quiet Sun have been taken to imply that some process 
preferentially enhances the helium resonance line intensities 
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with respect to other lines formed in the transition region 
in the quiet Sun, and that the enhancement is reduced in 
coronal holes. 

Current evidence suggests that in the quiet Sun the 
He II, and to a lesser extent the He I, resonance lines are 
formed mainly by collisional excitation (see Smith & Jor- 
dan 2002 - hereafter Paper I - for supporting references). 
Because the helium resonance lines have unusually large val- 
ues of W/kTe, where W is the excitation energy and Te is 
the electron temperature, their collisional contribution func- 
tions are sensitive to excitation by suprathermal electrons. 
He I and He II also have long ionization times compared with 
other transition region species, which may allow departures 
from ionization equilibrium in helium. Any process exposing 
helium ions to larger populations of suprathermal electrons 
than in equilibrium will tend to increase the collisional ex- 
citation rates of the helium lines, while lines with smaller 
W/kTe will be relatively unaffected. MJ99 reported that at- 
tempts to model such effects had been generally unsuccessful 
in explaining the helium resonance line emission relative to 
other transition region (TR) lines. While a conclusive expla- 
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nation is still lacking, some recent work appears promising 
(see Paper I). 

Fontenla, Avrett & Loeser (2002) have extended their 
earlier work (Fontenla et al. 1993), performing radiative 
transfer calculations for hydrogen and helium including the 
effects of mass-conserving flows as well as ambipolar diffu- 
sion and departures from ionization equilibrium. His prelim- 
inary results show improved agreement with observations, 
outflows (of 5-10 km s" 1 at T = 10 5 K) producing increases 
in intensity in the He I and He II resonance lines of up to 
an order of magnitude. Inflow models produce self-reversals 
in the He I 584. 3- A line, whereas observations show the line 
is more likely to be self-reversed in cell interior regions (e.g. 
MJ99) where the line shows a blue shift (Peter 1999). This 
implies that the models do not yet include all relevant pro- 
cesses. As the models do not include the effects of diffusion 
on any elements heavier than helium, it is also uncertain as 
to how successful they will be in explaining the helium line 
intensities relative to other transition region lines. 

Jordan (1975,1980) suggested two possible enhance- 
ment mechanisms dependent on the magnitude of the tem- 
perature gradient that might explain the coronal hole/quiet 
Sun contrast, Munro & Withbroe (1972) having found 
dT/dh to be an order of magnitude smaller inside coronal 
holes. The first involves the transport of helium ions up the 
steep TR temperature gradient, allowing them to be excited 
at electron temperatures higher than those in equilibrium. 
Jordan (1980) found that intensity enhancements of up to 
a factor of 5 could be produced in the He n 303. 8- A line 
in the quiet Sun, and a more detailed study by Andretta et 
al. (2000) suggests similar factors. A further investigation of 
the process using an improved treatment of the excitation 
times of the helium ions is reported in Paper I. The analysis 
was extended to include estimates of the effects on the He I 
584. 3- A and 537. 0-A lines. We conclude, as did Andretta et 
al. (2000) for the He n line, that current observations are 
consistent with the process accounting for at least part, and 
perhaps all, of the enhancement apparently required in the 
helium resonance lines in the quiet Sun. 

A complementary study of the other of Jordan's (1980) 
suggestions, enhanced collisional excitation by suprathermal 
electrons of non-local origin, is presented in this paper. Ra- 
diative transfer codes usually assume that the electron veloc- 
ity distribution functions (EVDFs) are Maxwellian through- 
out the model atmospheres, but observations of the solar 
wind (e.g. Scudder 1994) and numerical studies of the form 
of the EVDF in the solar transition region suggest that sig- 
nificant departures from Maxwellian distributions may oc- 
cur. This may result from some acceleration process main- 
taining a non- Maxwellian distribution in the low TR or chro- 
mosphere (Scudder 1992a; Vinas, Wong & Klimas 2000) or 
from the transport of high energy electrons from the upper 
TR and corona (Shoub 1983; Ljepojevic & Burgess 1990). In 
the second case at least, the influence of the suprathermal 
electrons would depend on the temperature gradient. 

The presence of an enhanced suprathermal electron 
population would simulate an increase in electron tem- 
perature in transitions with excitation energies well above 
the local thermal energy, and hence lead to larger excita- 
tion rates. Transitions with W smaller than the energy at 
which the suprathermal electron population becomes im- 
portant would not be affected significantly. Postulated non- 



Maxwellian EVDFs generally show departures from the local 
Maxwellian distribution which increase monotonically with 
energy, so that their effects on ionization rates can be even 
larger than on line excitation. Essentially the same effect 
limits the effectiveness of non-thermal transport of helium 
at temperatures much higher than those of normal line for- 
mation. In the case of collisional ionization by non-local 
suprathermal electrons the effect is a shift of the ioniza- 
tion equilibrium; increased collisional ionization rates raise 
the population of each ion at lower temperatures. Emission 
could therefore be increased because it would occur over a 
larger, more dense region. 

In a study of Si in emission line ratios, Pinfield et al. 
(1999) claim to have observed both an enhancement of the 
intensity of the line with the largest value of W/ kT and a re- 
duction in peak line formation temperature in the quiet Sun 
with respect to a coronal hole. In light of the above discus- 
sion, this is evidence of excitation by an enhanced suprather- 
mal population of electrons whose importance depends on 
the temperature gradient in the transition region. 

Although the same mechanism has often been suggested 
to explain the anomalous intensities of the helium resonance 
lines, few detailed calculations of the possible effects on he- 
lium have been made, perhaps because of the difficulty of 
solving the Boltzmann equation for the EVDF and the equa- 
tion of radiative transfer simultaneously. Shoub (1983) com- 
puted non-Maxwellian collision rates for the ionization of 
He I and He II and the excitation of the 584. 3-A and 303.8- 
A lines, but did not calculate line intensities for compari- 
son with observations. Anderson, Raymond & Ballegooijen 
(1996) computed intensities for the He II 303. 8-A line us- 
ing non-Maxwellian EVDFs, but did not include radiative 
transfer, and made no study of He I. 

In this paper I present the results of radiative trans- 
fer calculations for cases of both Maxwellian and non- 
Maxwellian collisional excitation. The former illustrate the 
problems found in attempting to model the formation of the 
helium resonance lines; the latter are used in a study of the 
observable effects that non-Maxwellian EVDFs might have 
on He I and He II lines. The calculations allow an extensive 
comparison of the predictions of the enhanced excitation 
mechanism with observations. The atomic and atmospheric 
models used in the radiative transfer calculations are de- 
scribed in Section 2, and the results of radiative transfer 
calculations performed with Maxwellian collision rates are 
presented in Section 3. Section 4 describes how the effects of 
non-local electrons are simulated. The results of these latter 
calculations are presented and discussed in Section 5, and 
the conclusions are summarized in Section 6. 



2 ATOMIC AND ATMOSPHERIC MODELS 

The radiative transfer modelling of the solar helium spec- 
trum was carried out using version 2.2 of the multi code 
(Scharmer & Carlsson 1985; Carlsson 1986). This code may 
be used to solve non-LTE problems in semi-infinite plane- 
parallel one dimensional model atmospheres. In order to 
calculate an emergent spectrum, the code requires a model 
atom, a model atmosphere, and a set of abundances (which 
are used to calculate background opacities - these are gen- 
erally computed assuming LTE, but in the calculations re- 
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ported here, the opacity due to hydrogen was computed in 
non-LTE). It also allows the specification of the coronal ra- 
diation field as a boundary condition at the upper edge of 
the model atmosphere. 

2.1 The model atom 

The atomic model used in this study comprises 29 bound 
levels of He I (lsnl ' L terms with principal quantum num- 
ber n < 5), six bound levels of He II with n < 6, and the 
He ill ground state. The atomic data used before modifi- 
cations were made to incorporate collisional excitation and 
ionization by non-Maxwellian EVDFs are listed below. 

The He I part of the model was initially based on the 30 
level He I model of Andretta & Jones (1997), but using newer 
and/or more complete data for some of the parameters, and 
was extended to include He II. Greater detail was included 
in the neutral stage, in order to facilitate investigation of the 
relationships between the 584. 3-A and 537. 0-A lines, which 
have been observed extensively with SOHO. Previous work 
by Hearn (1969) and Andretta & Jones (1997) suggested 
that the formation of these lines depends on a combination 
of collisional and radiative transitions between many lev- 
els, and is not necessarily dominated by either direct colli- 
sional excitation or photoionization- recombination (PR). It 
was therefore important to include as many levels and pro- 
cesses as possible in this part of the model. In He II most 
attention was paid to the 303. 8-A resonance line, which also 
appears in SOHO CDS spectra. 

The energies of the He I levels are taken from Martin 
(1987), averaging over the fine structure in the triplet lev- 
els. The fine structure in the triplet terms is not included, 
which is a valid approximation in the rate equations since 
the energy separation of the sub-levels is small enough for 
collisional transitions between them to keep their popula- 
tions in the proportions of their statistical weights. As noted 
by Andretta & Jones (1997), although this assumption may 
not be valid in the 2p 3 P° term, where the separation of the 
J-states is greatest, radiative transitions will also tend to 
populate the sub-levels in the same proportions in the con- 
ditions encountered in the atmospheric models used. The 
energies of the He II levels are from Kisielius, Berrington & 
Norrington (1996) for n < 4, from Bashkin & Stoner (1975) 
for n > 5. States of different I were summed over according 
to their statistical weights. 

Oscillator strengths for most of the allowed transitions 
in He I are provided by Drake (1996), while values for 
the few transitions not covered there are taken from Theo- 
dosiou (1987). The rate for the spin- forbidden electric dipole 
\s2p 3 Pi - Is 2 1 S transition is taken from Drake & Dalgarno 
(1969), and the two photon rate for the ls2s 1 S - Is 2 1 S 
transition is also included, from Bassani & Vignale (1982). 
The He II oscillator strengths are taken to be hydrogenic, 
using the results of Wiese, Smith & Glennon (1966). 

Stark broadening parameters for He I are from Dimitri- 
jevic & Sahal-Brechot (1984,1990) where available, with the 
remainder being estimated using the formula of Freuden- 
stein & Cooper (1978). The Stark widths of Griem (1974) 
are used for He II. Van der Waals broadening is treated us- 
ing the tables of Deridder & Van Rensbergen (1976), taking 
advantage of alterations made to multi by Rowc (1996). 

For most of the levels of He I photoionization cross- 



sections are drawn from the Opacity Project database 
(Seaton 1987); they were calculated by Fernley, Taylor & 
Seaton (1987). The broadest, most prominent resonances 
in the cross-sections are represented approximately in the 
model, particularly in the resonance continuum, as there is 
some variation of the coronal radiation field in the wave- 
length range of the resonances. The 4 3 P, 4 1,3 F, 5 3 P, 

5 1,3 F, and 5 1,3 G levels of He I are treated as hydrogenic, 
as are the photoionization cross-sections for He II (Menzel 

6 Pekeris 1935). 

Photoionization is important in establishing the ioniza- 
tion balance between all three stages of helium and is im- 
portant in the formation of some of its lines. The effects of 
coronal radiation shortwards of 504 A are therefore included 
in the calculation of ionization and recombination rates. The 
coronal illumination was represented using the solar EUV 
irradiance model of Tobiska (1991,1993). The spectrum of 
radiation at minimum coronal activity was used to give in- 
coming intensities at the solar surface representative of the 
average quiet corona. Photoionization of He I occurs directly 
from the ground state, by coronal radiation, and by a two- 
stage process, in which He I triplet states, collisionally ex- 
cited from the ground, are rapidly photoionized (directly 
or by successive photoexcitations followed by ionization) by 
the photospheric radiation field. The importance of this two- 
stage process was recognized in early calculations by Athay 
& Johnson (1960) and Hearn (1969), and by Andretta & 
Jones (1997), who emphasized how the singlet and triplet 
systems of He I have quite different roles in the ionization 
balance. Collisional excitation of singlet levels tends to lead 
more often to decays back to the ground, to which the sin- 
glet levels are strongly connected by the allowed radiative 
transitions lsnp 1 P° - Is 2 1 S. In the triplet system, the 
metastable 2 3 S level acts like a ground state, and recombi- 
nations to higher triplet levels tend not to lead to cascade 
decays, but more often result in re-ionization by the pho- 
tospheric radiation field. Recombinations to singlet states 
generally result in radiative cascades to the ground. 

Collisonal excitation rates between bound levels of He I 
are calculated, where possible, using the collision strengths 
of Lanzafame et al. (1993), as these are valid over a large 
temperature range. Collision rates in the resonance series are 
among those calculated using these data. Collision strengths 
for most of the remaining transitions in He I are provided 
by Sawey & Berrington (1993), and gaps are filled using 
data from Mihalas & Stone (1968) and Benson & Kulan- 
der (1972). Some of these early data are expected only to 
be correct to an order of magnitude. Collision strengths for 
bound-bound transitions in He II are drawn from Aggar- 
wal et al. (1992) for n < 3 and from Aggarwal, Berrington 
& Pathak (1991) for transitions involving n — 4. The re- 
maining He n collision strengths are calculated using the 
expression given by Mihalas & Stone (1968). 

Collisional ionization rates from all levels of He I and 
He II are calculated using expressions from Mihalas & Stone 
(1968). Older data have been used in preference to newer 
information because the more recent work is not general 
enough. Bray et al. (1993), for example, give the ionization 
cross-section for only the ground state of He II; the equiva- 
lent cross-section provided by Mihalas & Stone (1968) is in 
fact very similar, but they also give expressions for ioniza- 
tion from excited states. 
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Figure 1. Temperature versus column mass density (m) in the 
VAL C (solid) and D (dashed) model atmospheres. 

Dielectronic recombination to He I, for both Maxwellian 
and non-Maxwellian electron distributions, was investigated 
in detail (Smith 2000), but was found to have very little 
effect on the formation of the 584. 3-A and 537.0-A lines. 



2.2 Model atmospheres 

Since the detailed physical processes responsible for heat- 
ing the chromosphere, transition region, and corona remain 
unknown, most of the atmospheric models suitable for use 
in radiative transfer calculations are semi-empirical. Their 
temperature structures are chosen to reproduce observations 
of spectral lines and continua, rather than being derived 
from the energy balance ab initio. The VAL model atmo- 
spheres (Vernazza, Avrett & Loeser 1981, hereafter VAL), 
are successful in many respects and are the most often used 
solar models of this type. The models were constructed to 
fit six brightness components (A-F) observed by the Har- 
vard instruments on Skylab. The C component represents 
the average quiet Sun, and the VAL C model was initially 
adopted in the present work. The VAL D model, represent- 
ing the average quiet Sun network, was also tested against 
observations. The run of temperature against column mass 
density in both models is shown in Figure 1. 

Fontenla et al. (1993) modified the VAL model atmo- 
spheres in an effort to improve the physical basis of the mod- 
els. These new (FAL) models have semi-empirical chromo- 
spheric structure based on the VAL models, but have tran- 
sition regions derived from energy balance considerations 
including the the effects of ambipolar diffusion of H and He. 
These models produce results consistent with observations 
of the hydrogen lines without need for the plateau at about 
2 x 10 4 K in the VAL models (see Figure 1), but do not re- 
produce observed helium resonance line intensities well. The 
FAL models neglect the effects of turbulence (except in the 
equation of hydrostatic equilibrium), and it is not clear to 
what extent the steady state diffusion flows in the models 
are present in the dynamic solar atmosphere. The models 
also do not include the effects of diffusion on any elements 
heavier than helium, lines of which must be compared with 
the helium lines in any explanation of the formation of the 
latter. These factors suggest that the representation of the 



transition region in the FAL models is not necessarily su- 
perior to that in the semi-empirical VAL models. For these 
reasons, the VAL models were used in preference to FAL in 
the present work. 

There are some problems with the VAL models. Calcu- 
lations using the VAL C and D model atmospheres produce 
intensities in the C II lines around 1335 A about a factor of 
six larger than the values observed by MJ99 (Smith 2000). 
This appears to be due to the influence of the 'Lyman a 
plateau' in the VAL models at about 2 x 10 4 K, which was 
included in order to put enough material in the appropriate 
temperature range to reproduce the observed intensities of 
lines of H and He, particularly the very important H Lye* 
line. The plateau produces a local peak in the EMD in the 
temperature range of the formation of both the He I reso- 
nance series and the C II lines. MJ99 found that this peak 
also causes Si ill line intensities to be over-produced by a 
factor of at least 2 with respect to their observations. In 
view of this problem with the VAL models, and in order to 
make comparisons with MJ99 more meaningful, new model 
atmospheres were constructed to be consistent with emis- 
sion measures derived from observed intensities of TR lines 
of elements other than helium. 

The new models were chosen to represent the network 
rather than cell interior regions. Models of the magnetic field 
generally show flux tubes confined to the network bound- 
aries at photospheric levels expanding in the chromosphere 
and TR to fill the atmosphere in the corona (see e.g. Gabriel 
1976). Such models suggest that network regions may be 
more reasonably modelled with radial magnetic fields than 
cell interiors. The effects of non-local electrons would be 
expected to be greatest where the magnetic field extends di- 
rectly into the corona. Also, the VAL chromospheric model 
may not be valid in cell interiors, as the time-dependent 
models of Carlsson & Stein (1995) suggest that a classical 
chromosphere with an outward-increasing temperature may 
not exist in non-magnetic internetwork regions. 

The multi code requires, as input for a (static) model 
atmosphere, the run of mass column density, temperature, 
electron density, and microturbulent velocity. Below T ~ 10 4 
K these figures were taken directly from the VAL D network 
model. The transition region and coronal parts of the mod- 
els were developed from emission measure distributions us- 
ing methods described by Philippides (1996) and McMurry 
(1997). The EMDs used to derive these parts of the new 
models are described in detail in Paper I and only brief de- 
tails are repeated here. In the upper TR (above log Te = 
5.3), the EMDs were derived assuming energy balance be- 
tween radiation losses and the divergence of the classical 
conductive flux from the corona. The coronal emission mea- 
sures were chosen to fit the EMD to MJ99's observations of 
Mg IX and Mg x. As discussed in Paper I, this requires pres- 
sures higher than were observed by MJ99 in the network, but 
similar to those in the VAL models. Below log Te = 5.3, the 
lower transition region of model S was derived directly from 
MJ99's network EMD. As explained in Paper I, the form 
of the EMD is uncertain, as the CDS and SUMER lines could 
not be observed at the same locations at the same time. A 
second model was therefore constructed. Model X was made 
using the EMD derived from HST observations of the G8 V 
star £ Boo A, scaled to the minimum of the MJ99 network 
EMD (see Paper I) . The lower transition region parts of the 
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Figure 2. Temperature versus column mass density in network 
models S (dashed) and X (dotted) in the region where they differ 
from the VAL D model (solid). 



two models were derived assuming hydrostatic equilibrium 
(including a contribution to the pressure from turbulent mo- 
tions) using methods discussed by Jordan & Brown (1981) 
and more recently by Harper (1992). 

The resulting transition region models were grafted sep- 
arately on to the VAL D model chromosphere. Running 
multi in hydrostatic equilibrium with a 9 level hydrogen 
atom produced new values of x(T) = N^/Ne and of the 
height at the base of the transition region, which were used 
to improve the transition region parts of the models. Self- 
consistent models were found by iterating this process to 
convergence in x(T). The resulting models are shown in Fig- 
ure 2, which shows how the temperature plateau in the VAL 
D model is suppressed in the new network models. The full 
models used in the calculations are given in an appendix in 
Tables 7 and 8. 

The elemental abundances used in the calculations were 
the photospheric abundances of Grevesse, Noels & Sauval 
(1992) and Anders & Grevesse (1989). The helium abun- 
dance is taken as Njn e /Njn = 0.098 and is assumed to be 
constant through the atmosphere. Calculations by Hansteen, 
Leer & Holzer (1997) predict that the helium abundance 
could vary significantly in the solar atmosphere. Their mod- 
els show a decreased helium abundance in the transition re- 
gion, and so cannot explain the apparent enhancement of 
the helium line intensities. 



3 RESULTS USING MAXWELLIAN 
COLLISION RATES 

In this Section the results of calculations assuming purely 
Maxwellian electron distributions are discussed, and com- 
pared with the observations of MJ99. Calculations were per- 
formed both with and without the quiet coronal illumination 
described in Section 2 in order to investigate the importance 
of photoionization-recombination (PR) in the formation of 
the helium resonance lines. The discussion focusses on the 
He n resonance line at 303.8 A and the first two lines of 
the He I resonance series, at 584.3 A and 537.0 A. The He I 
triplet lines (e.g. 10830 A) are not discussed here; they were 



Table 1. Quiet Sun intensities and line ratios observed by 
Macpherson and Jordan (1999). For averages over the network 
see also Jordan et al. (2001). 



Region: 


Strong 


Typical 


Cell 


Average 




network 


network 


interior 


quiet Sun 


A 




Observed intensity 




(A) 




(erg cm 2 


s 1 sr 1 ) 




584.3 


613 


613 


380 


487 




±93.0 


±114 


±126 




537.0 


72.8 


70.9 


44.7 


56.9 




±10.5 


±12.3 


±17.1 




303.8 


10050 


8289 


5066 


6654 




±1249 


±1216 


±1540 




/(A)//(584.3 A) 


537.0 


0.119 


0.116 


0.118 


0.117 


303.8 


16.4 


13.5 


13.3 


13.7 



discussed in detail by Andretta & Jones (1997), and nothing 
has been found that contradicts their analysis. 

The intensities observed by MJ99 in different regions of 
the quiet Sun are given in Table 1, with an 'average quiet 
Sun' intensity calculated for comparison with results from 
the VAL C model. On the basis of a comparison of the 
intensity classifications of MJ99 and VAL, this average is 
weighted as 0.54 x 'cell interior' intensity + 0.40 x 'typ- 
ical network' intensity ± 0.06 x 'strong boundary' inten- 
sity. The intensities given by MJ99 were derived using the 
calibration in the CDS software at the time of the observa- 
tions (1997), with the modifications suggested by Landi et 
al. (1997). Jordan et al. (2001) show how these intensities 
would be changed by adopting the calibration derived by 
Brckkc et al. (2000). When the most recent calibration is 
used, the main effect is a reduction of the He II 303. 8- A line 
intensities by a factor of 2.2. 

3.1 Quiet coronal illumination 

The intensities of the resonance lines computed with quiet 
coronal illumination are given in Table 2. As found by An- 
dretta & Jones (1997) in their study of the He I lines, ex- 
citation of the 2 1 P level occurs directly from the ground 
by collisions and, at least as rapidly, through the excitation 
of levels with strong allowed radiative transitions to 2 1 P. 
Over much of the region of line formation, the most impor- 
tant populating process is by photoexcitation from 2 1 S, 
which itself is mainly excited by collisions from the ground. 
Allowed transitions from other singlet states excited from 
the ground, particularly 3 S and 3 1 D, are also important; 
net rates for the two-step processes are comparable to direct 
collisional excitation of 2 1 P from the ground. The contribu- 
tion function peaks at log Te — 4.45. Collisional excitation 
(direct and indirect) dominates above this temperature. At 
lower temperatures (4.0 < log Tq < 4.4) collisional excita- 
tion is less important than radiative recombination cascades 
through the singlet levels (3 1 S, 3 1 D etc.). 

Using the VAL C model atmosphere, the calculated in- 
tensity of the 584. 3-A line is smaller (by about 20 per cent) 
than the quiet Sun average of MJ99's observed intensities 
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Figure 3. He I 584. 3-A line profiles computed using the VAL C 
(top) and network S (bottom) models with quiet coronal illumi- 
nation. 



(see Tables 1 and 2) . Output from the code suggests that the 
line is optically thick, with an optical depth at line centre be- 
tween 10 and 100, but it is effectively thin, as found by Hearn 
(1969); i.e. almost all photons created in the line escape the 
atmosphere. The computed line profile has a deep central 
self-reversal (see Figure 3) . At the spectral resolution of the 
SUMER instrument (about 0.04 A; Wilhelm et al. 1995), such 
a reversal would be obscured by instrumental broadening. A 
convolution of the computed profile with a Gaussian instru- 
mental width of 0.04 A shows a much smaller self-reversal 
(see Figure 4) . This convolved profile is more consistent with 
observations of the quiet Sun. Using SUMER, MJ99 found 
evidence for a small self-reversal in cell interior regions, but 
none in the network. Peter (1999) found the mean quiet Sun 
profile to be flat-topped. In coronal holes, where the line has 
a larger optical depth, both Peter (1999) and Jordan et al. 
(2001) observed self-reversed profiles. The computed width 
of the line (FWHM ~ 0.13 A) is similar to those obtained 
in calculations by Andretta & Jones (1997), Fontenla et al. 
(1993; 2002). It also compares well with the mean width of 
0.14 A observed by Doschek, Behring & Feldman (1974), 
and the widths < 0.13 A observed by Peter (1999). 

The 537. 0-A line forms by similar processes to the 584.3- 
A line, with contributions to the population of the upper 
level, 3 1 P, from direct collisional excitation and radiative 
transitions from other levels populated by a combination of 



collisional excitation and recombination. The most impor- 
tant intermediate levels in this case are 2 1 S, 4 1 S and 4 1 D. 
Unlike in the 584. 3- A line, direct recombination to 3 1 P 
is more important than cascade recombination. The overall 
contribution function is similar to that of the 584.3 A line, 
but peaks at slightly lower temperatures (log Te ~ 4.40). 
The 537.0-A line has an optical thickness smaller than than 
that of the 584. 3-A line, but is less likely to be effectively 
thin, as explained by Hearn (1969). The optical depth at line 
centre is of order 10 in calculations with the VAL C model 
atmosphere. The computed intensity of the line is about half 
the quiet Sun mean of MJ99's observations. 

The formation of the He II resonance line is domi- 
nated by direct collisional excitation, its contribution func- 
tion peaking at a temperature of log Te = 4.9. There is, 
however, a small but not negligible contribution to the pop- 
ulation of the n = 2 level from radiative cascades from higher 
levels. The levels with n > 2 are similarly populated by col- 
lisional excitation at high temperatures and recombination 
cascades at low temperatures. Recombination becomes in- 
creasingly important as n increases (see below for a brief 
discussion of the 1640.4 A multiplet). The 303.8 A line is 
marginally optically thick, with optical depth at line centre 
greater than 1 but less than 10. This agrees with separate 
approximate calculations, using the parameters of the radia- 
tive transfer models, predicting optical depths in the He II 
line of order 5 (and of order 100 in the He I 584. 3-A line) . 
The line is certainly effectively optically thin. In results from 
the VAL C model, the computed line profile is approximately 
Gaussian, and shows no sign of self-reversal (see Figure 7 in 
Section 5). The line width (FWHM ~ 0.08 A) is somewhat 
smaller than the value of 0.10 A observed by Doschek et al. 
(1974) and Cushman & Rense (1978), but is very similar 
to that computed by Fontenla et al. (1993). Fontenla et el. 
(2002) predict greater widths if significant flows are present. 
The integrated intensity calculated for the VAL C model is 
a factor of more than 3 smaller than the quiet Sun average 
of MJ99's observations (or only a factor of 1.5 smaller when 
the most recent CDS calibration is used), but is about twice 
that computed by Fontenla et al. (1993). 

MJ99 suggested that the appearance of the network in 
the helium resonance lines (broader and with lower contrast 
than in other lines formed in the low TR) could be explained 
by their large optical depths. If the network structures are 
thicker in the line of sight than perpendicular to it, pho- 
tons would more likely escape from the edges of the network 
into cell interiors. If then resonantly scattered back into the 
line of sight, this would apparently broaden the network, 
and might explain why the helium lines appear to be most 
enhanced with respect to other TR lines in the cell interi- 
ors. Both the He I and He II resonance lines were found to 
be optically thick in the calculations described here, but as 
these calculations were performed in the plane-parallel ap- 
proximation, detailed conclusions about scattering from the 
network cannot yet be drawn, but await radiative transfer 
calculations with two-component atmospheric models. 

Table 2 also gives the intensities of the helium lines 
computed using the VAL D network model and the new 
network models S and X. The 584. 3-A line intensity com- 
puted using the VAL D model is close to the observed mean 
network value, but the ratio 7(537.0 A)/7(584.3 A) is again 
smaller than observed. The computed 303. 8-A line intensity 
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Table 2. Integrated intensities computed for the VAL C and D 
and network S and X models with quiet coronal illumination. 



Model atmosphere: 


VAL C VAL D S 


X 


A 




Computed intensity 




(A) 




(erg cm~ 2 s _1 sr _1 ) 




584.3 


406 


624 146 


81.2 


537.0 


31.2 


50.7 17.1 


6.78 


303.8 


2020 


6400 1163 


988 


/(A)//(584.3 A) 


537.0 


0.077 


0.081 0.117 


0.083 


303.8 


4.98 


10.3 7.97 


12.2 



1.5x 10 



1.0x10 



(5 5.0x10 : 




583.6 583.8 584.0 584.2 584.4 584.6 584.8 585.0 

Wavelength (A) 



is only 25 per cent smaller than was observed by MJ99, but 
adopting the most recent calibration, the computed inten- 
sity is larger than observed by a factor of 1.7. The VAL 
D model also predicts other low TR lines (e.g. of C n) to 
be stronger than observed by significant factors. Models S 
and X were constructed to be more consistent with the ob- 
served intensities of TR lines other than those of helium. 
Significant changes in the lower atmospheric structure with 
respect to the VAL models lead to changes in the relative 
importance of the various processes responsible for the for- 
mation of the He I lines. Although the electron pressures are 
higher in the new network models than in the VAL C model, 
the intensities of the He I lines are smaller owing to the ab- 
sence of the temperature plateau which is present in the 
VAL model. This reduces the amount of material present in 
the new models in the temperature range where collisional 
excitation (direct and indirect) is effective, which is limited 
from below by the temperature dependence of the excita- 
tion rates and from above by the ionization of He I. Optical 
depth at the 584. 3-A line centre is still greater than 10. PR 
provides a relatively greater contribution to the total excita- 
tion rate than in the VAL calculations, particularly in model 
X. 

The intensities computed for the 584. 3-A and 537. 0-A 
lines using model S are a factor of about 4 smaller than 
observed network intensities, but the ratio of the computed 
intensities is close to the mean quiet Sun network value of 
0.1 16 ±0.015 found by MJ99. That value assumed the Landi 
ct al. (1997) calibration of CDS; using instead the most recent 
calibration reduces the ratio to 0.105 ± 0.014. While this 
is matched less well by the computed ratio, model S still 
produces a better fit than the other model atmospheres. 
The computed profile of the 584. 3-A line is also quite close 
to that observed, showing only a small self-reversal that is 
completely obscured by instrumental broadening of 0.04 A 
(see Figures 3 and 4). The intensities computed using model 
X are smaller than observed by factors of 7.5 for the 584. 3-A 
line and 10 for the 537. 0-A line, owing to the lower emission 
measure derived from the scaled £ Boo A distribution used 
in the lower TR. The ratio of the line intensities is smaller 
than observed. 

The value of the intensity ratio 7(537.0 A)/7(584.3 A) 
found using the VAL C model atmosphere is very similar to 
that obtained by Andretta & Jones (1997) in an equivalent 
calculation; they found the ratio to be larger in models in 
which the VAL C plateau was removed. Similarly, in mod- 
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Figure 4. He I 584. 3-A line profiles computed using the VAL 
C (top) and network S (bottom) models (shown in Figure 3), 
convolved with a Gaussian function representative of the instru- 
mental broadening of the SUMER instrument. 



els S and X, where the plateau is suppressed, the value of 
the computed ratio is larger than in results from the VAL 
models. The reduction in the size of the line forming region 
lowers the optical depth in the 537.0-A line (The 584.3-A 
line is effectively thin in all the atmospheric models tested 
here, but the 537.0-A line is not effectively thin in the VAL 
C calculations). Thus, although fewer photons are created in 
the 537.0-A line in the network models, a greater proportion 
escape, allowing the intensity ratio 7(537.0 A)/7(584.3 A) to 
increase compared with the VAL C case. The ratio is smaller 
in results from model X than in those from model S, al- 
though the region where collisional excitation of the lines is 
important is narrower in model X than in model S. In model 
X, the line formation has a relatively greater contribution 
from PR from the thicker region below T e ~ 1.8 x 10 4 K, 
where the optical depth in the 537.0-A line is greater. In re- 
combination cascades, levels linked to both 2 1 P and 3 X P by 
allowed radiative transitions generally decay to the former 
more rapidly. 

In the region of formation of the He II resonance line the 
new models have electron densities smaller than in the VAL 
D model, but have steeper temperature gradients, reducing 
the extent of the emitting region and thus the emergent 
intensity. The line still has optical depth greater than 1. The 
line intensities computed using models S and X are factors 
of 7 and 8 smaller than the mean observed network value 
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respectively. These factors are reduced to 3-4 when the most 
recent CDS calibration is assumed. The temperature gradient 
is slightly larger and the density smaller in model X than 
in model S, explaining the lower emission in model X. The 
ratio 1(303.8 A)/J(584.3 A) found using model X (12.2) is 
the closest of those computed to the observed network ratio 
of 13.5. using the most recent CDS calibration, the observed 
ratio is about a factor of 2 smaller, and the ratio found 
using model S is a closer match. The computed profiles of 
the 303.8 A line, with FWHM ~ 0.1 A, are broader than in 
the VAL C results because of higher turbulent velocities in 
the regions of line formation in the new models. 

The separate components of the He II Balmer a mul- 
tiplct at 1640.4 A are not included explicitly in the cur- 
rent atomic model, but the total predicted intensity may 
be compared with observations. In calculations using the 
VAL C model with quiet coronal illumination recombination 
is more important than collisional excitation in the forma- 
tion of the multiplet, and the computed intensity is almost 
a factor of 6 smaller than the value of 88 erg cm~ 2 s _1 
sr _1 observed in the quiet Sun by Kohl (1977). Models S 
and X produce even poorer matches to observations, pre- 
dicting intensities at least an order of magnitude too small. 
Wahlstr0m & Carlsson (1994) found a slightly smaller dis- 
crepancy than in the results here, of a factor of about 4, but 
this is still a significant problem that requires explanation. 

3.2 Variations in coronal illumination 

Andretta & Jones (1997) investigated the effects on the He I 
lines of varying the intensity of the coronal radiation field by 
multiplying the assumed 'quiet' illumination by a constant 
factor at all wavelengths. They found that both the 584.3- 
A and 537. 0-A lines increased in intensity with increasing 
coronal illumination owing to increased contributions to line 
formation from PR. The results found here using the VAL 
C model agree with those of Andretta & Jones (1997), and 
calculations using models S and X show trends similar to 
results from their 'no plateau' models. The effects on He I of 
increasing the coronal intensity are therefore not discussed 
in detail (but are given in Smith 2000). Andretta & Jones 
(1997) did not, however, study lines of He II. 

Somewhat surprisingly, the He n 303. 8-A line, whose 
formation in all cases is dominated by collisional excitation, 
generally shows a decrease in intensity as the coronal radi- 
ation is increased. This appears to be due to a shift of the 
peak He n ionization fraction to lower temperatures, reduc- 
ing the collisional excitation rate relative to the quiet corona 
case. This outweighs any increase in the small contribution 
from recombination to the line, although when the coronal 
irradiance is reduced to zero in model X the computed in- 
tensity of the He II line is also slightly reduced. Given that 
spatial variations of the intensities of the He I and He II 
resonance lines are very similar (e.g. MJ99), and certainly 
not anti-correlated, this result supports the contention that 
PR by coronal radiation does not dominate the formation 
of both lines. 

The model atmospheres considered in this work are not 
appropriate to coronal hole regions, where the density is 
a factor of 2 or more smaller than in the quiet Sun, and 
the temperature gradient is much shallower. Removing the 
coronal illumination included in the radiative transfer cal- 



culations does however indicate how the lower photoionizing 
flux in coronal holes might affect the helium resonance lines. 
The intensities of the lines are observed to be a factor of 1.5 
- 2 smaller in coronal holes than in the quiet Sun (Peter 
1999; Jordan et al. 2001). Much smaller reductions in the 
intensities of the He I lines are computed when the coronal 
radiation incident on the top of the VAL C model is removed. 
In results from models S and X, where PR is more important 
in forming the He I lines, more significant effects are seen. 
The computed intensity of the 584.3 A line in model S is re- 
duced by a factor of 1.3 when the quiet coronal illumination 
is removed, while with model X this causes the line inten- 
sity to fall by a factor of almost 4. The observed increase 
in the value of the He I line ratio 1(537.0 A)/J(584.3 A) in 
coronal holes is understandable in terms of the reduction of 
coronal illumination. Such a trend is seen in the results from 
all of the models, and the value of 0.134 found from model 
S with zero coronal illumination agrees quite well with the 
mean value of 0.130 ± 0.014 found in coronal hole observa- 
tions by Jordan et al. (2001), although this observed value 
is decreased to 0.117 when the most recent CDS calibration 
is adopted. 

In none of the models tested here is the He n resonance 
line reduced significantly in intensity by the removal of coro- 
nal illumination. As this line is formed mainly by collisional 
excitation its intensity would be expected to depend to a 
greater extent on changes in the density and temperature 
gradient in coronal holes than on changes in the coronal il- 
lumination. Given that the He I and He n resonance lines 
respond differently to significant changes in coronal illumi- 
nation, the observation that the ratio of the lines changes 
little between coronal holes and the quiet Sun (Jordan et al. 
2001) again implies that this is not the factor controlling the 
changes in absolute intensity. 



4 NON-MAXWELLIAN COLLISION RATES 

The calculations reported above show that none of the model 
atmospheres tested gives a good match to all aspects of 
the observations. Collisional excitation by non-Maxwellian 
EVDFs was investigated as a possible explanation of these 
discrepancies. Modifications were made to the radiative 
transfer code to simulate the excitation and ionization of 
the helium atom by EVDFs with enhanced suprathermal 
tails of a form that has been suggested to exist in the solar 
transition region. 

4.1 Transition region electron distributions 

In order to calculate the effects on collisional excitation 
and ionization rates of a given non-Maxwellian EVDF, one 
needs to know its shape as a function of temperature. 
Strictly, an electron temperature cannot be defined for a 
non-Maxwellian EVDF. In the distributions postulated to 
exist in the transition region, however, the bulk of the elec- 
trons have a nearly Maxwellian velocity spectrum, which is 
used to define the temperature. 

In any plasma, the EVDF is a solution to the Boltzmann 
(or Fokker-Planck) equation: 

dt dx rn'dv \dt) coll 
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where / is the EVDF, v and x are electron velocity and 
position, F is the total force acting on the electron, and the 
term on the right hand side is the rate of change of / with 
time due to collisional redistribution of electrons in veloc- 
ity space. The Maxwellian distribution is the homogeneous, 
steady state solution found in a gas in thermodynamic equi- 
librium, which is independent of x and t. 

Shoub (1983) showed that, although the solar TR is 
weakly inhomogeneous (the mean free path of thermal elec- 
trons is at all points small relative to the temperature and 
density gradient scale lengths), the Spitzer-Harm (1953) so- 
lution is invalid at high velocities, where the mean free path 
of an electron increases as the fourth power of its velocity. 
Studies of the transition region EVDF in the case where the 
Knudsen parameter (the ratio of the electron mean free path 
to the scale length) is large have been made in one of two 
approximations. 

Some more recent work has focused on the process of 
'velocity filtration,' proposed by Scudder (1992a,b). This 
postulates a non-Maxwellian EVDF with an over-population 
of suprathermal electrons at the base of the transition re- 
gion, which have large enough velocities to climb the gravi- 
tational gradient into the corona, there forming the bulk of 
the plasma at T ~ 10 K. Studies of this process assume that 
the Knudsen parameter is of order one or greater, so that 
the electron fluid is effectively collisionless, and the RHS of 
equation (1) is zero. Anderson et al. (1996) found that col- 
lisionless velocity filtration is inconsistent with observations 
of transition region lines, and suggested that the collision 
term in the Boltzmann equation should not be neglected. 

In a different context, others (e.g Shoub 1983) have con- 
sidered a high velocity form of the Boltzmann equation in- 
cluding a collision term due to Landau (1936) describing the 
interaction of particles under inverse-square Coulomb forces. 
Solutions are computed for a model transition region with 
a prescribed temperature profile derived semi-empirically 
from observations or energy balance arguments. Such calcu- 
lations have been performed both for plane-parallel model 
TRs (Shoub 1983; Ljepojevic & Burgess 1990) and for coro- 
nal loop models (Ljepojevic & MacNeice 1988,1989). The 
heating of the corona is assumed to occur by some unspec- 
ified mechanism, and the coronal temperature is used as a 
boundary condition. 

Whereas in the velocity filtration models a non- 
Maxwellian EVDF is assumed to exist as a boundary or 
initial condition, in the collisional models non-Maxwellian 
EVDFs are not assumed to exist a priori, but their existence 
and form in the transition region are derived as results of 
the steep temperature gradient. Given that the existence of a 
steep temperature gradient in the TR is confirmed by many 
observations, whereas independent processes by which non- 
Maxwellian EVDFs may be maintained are at present only 
postulated, the collisional calculations seem to rest on more 
solid assumptions. For this reason, and because Anderson 
et al. (1996) cast doubt on the collisionless approach, the 
present work is based on calculations of the EVDF using 
the Landau equation. 

Numerical solutions of the Landau equation in the TR 
result in EVDFs, averaged over pitch angle, resembling (very 
slightly underpopulated) local Maxwellians at low veloci- 
ties, with more heavily populated tails diverging from the 
Maxwellian at a few times the thermal velocity. The high 



velocity tail of the distribution at a point in the lower TR is 
almost wholly populated by electrons originating from the 
near-thermal parts of distributions present at higher tem- 
peratures. The temperature gradient provides an excess of 
high energy electrons moving downwards; as their speed- 
dependent mean free path varies as v 4 , these electrons arc 
more influential than the similar excess of low energy elec- 
trons moving up the gradient. 

Shoub (1982) obtained analytical solutions to the Boltz- 
mann equation using a linearized form of the Bhatnagar- 
Gross-Krook (1954) (BGK) model for the collision term. 
Using this method, he found that the angle-averaged dis- 
tribution functions obtained were in reasonable agreement 
with those found in his numerical work, with very similar en- 
hanced suprathermal tails. The analytical solutions for the 
angle- averaged distribution in the lower TR can be accu- 
rately approximated by a power law over a wide range of 
suprathermal velocities. 

The distribution tested by Anderson et al. (1996) which 
most closely approximated the results of Shoub (1982,1983) 
was the k distribution (as suggested by Scudder 1992b) . An- 
derson (1994) suggested the use of a Maxwellian EVDF with 
a power law tail attached above some critical speed (as in- 
troducing a collision term makes the low speed part of the 
EVDF more nearly Maxwellian). This is actually a much 
better approximation to the collisional results of Shoub 
(1982,1983) and Ljepojevic & Burgess (1990), in which the 
transition from the Maxwellian form to a power law tail is 
generally sharper than in a n distribution. For these reasons 
I have chosen to use this Maxwellian plus power law EVDF 
in my calculations of enhanced collision rates. The calcula- 
tions are not formally self-consistent, but assume that the 
EVDF throughout the model atmospheres used in the ra- 
diative transfer calculations is of this form when calculating 
the collisional excitation and ionization rates for the model 
helium atom. 

4.2 Parametrization of the EVDF 

The way in which the EVDF is parameterized in the present 
work was influenced by the analytical BGK calculations 
made by Shoub (1982). I take the angle- averaged EVDF to 
be locally Maxwellian below a velocity £tail> w ith a power 
law decline at higher £, where £ is a dimensionless velocity 
v/ v th' an d v th ^ s t ' le thermal velocity ^/2fcT/m. A fea- 
ture of Shoub's (1982) solutions which is common to simi- 
lar work is that the value of £ marking the departure from 
near-Maxwellian is only very weakly dependent on the tem- 
perature defined by the Maxwellian bulk of the electrons, 
so as a first approximation I take £^ail to be constant with 
respect to T. 

The slope of the power law tail found by Shoub (1982) 
was determined by his choice of model atmosphere. This was 
an isobaric slab of thickness L in which energy is the trans- 
ferred by (classical) thermal conduction is constant. The cor- 
responding temperature profile takes the form 

T(z) 7/2 = T c 7/2 + (T^ 2 - Tc /2 )z/L (2) 

where Tc and are temperatures characterizing the in- 
coming EVDFs at the lower and upper boundaries of the 
atmosphere respectively. This choice of atmosphere is an 
obvious simplification. In the upper TR very little heating 
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is required to balance local radiation losses, and this can be 
provided easily by the net conductive flux. In the lower TR, 
however, this energy balance condition is not consistent with 
the rise in the emission measure distribution below T ~ 10 5 
K. Although the model atmosphere Shoub (1983) used in 
his numerical work allows for radiative power loss, the mod- 
els used in both papers represent the low TR poorly, hav- 
ing steeper temperature gradients than suggested by semi- 
empirical models. This failing may not have a major effect 
on the accuracy of the derived EVDFs in the lower TR, as 
the form of the suprathermal tail is more dependent on the 
temperature structure in the upper TR. In his analytical so- 
lutions, Shoub (1982) found that at T = 2 x 10 4 K, electrons 
with £ > 5 come mainly from T > 10 5 K, and at T = 8 x 10 4 
K, electrons with £ > 4 come from T > 2 x 10 5 K. 

Shoub's (1982) assumptions resulted in a power law de- 
cline in the tail of the EVDF varying as v~ :n/7 (or £~ 31/7 )- 
He generalized his results, however, to allow for a variation 
of the exponent r = 7/2 in equation (2). The corresponding 
index in the angle-averaged EVDF is — (12/r+l). This nota- 
tion has been used in the power law tails in the EVDFs used 
here. A smaller value of r, here designated ^ a ij, corresponds 
to a shallower temperature gradient low down and a steeper 
one high up, and to a steeper power law in the tail of the 
EVDF. As stated, Shoub's (1982) results have r tail = 3.5 
and are best represented in the formulation used here by tak- 
ing 3.5 < £ ta jj < 4.0. Both r^ a jj and £tail are taken to be free 
parameters in the calculations reported here, and ranges of 
values were investigated (see Figure 5). The parameter space 
was defined by the results of Shoub (1982,1983) and Ljepo- 
jevic & Burgess (1990). Calculations by the latter using two 
different model atmospheres produced EVDFs in the lower 
TR that may be approximated by taking 4.0 < Ctail — ^-0 
and 2.0 < r ta y < 4.0. Their solutions again showed £ t aa 
to be almost independent of temperature for a given model, 
except at low temperatures (T ~ 2 x 10 4 K). They found 
£tail ^° increase m this region, suggesting that the enhanced 
suprathermal tail may not persist right down to the bot- 
tom of the TR in models with more realistic (i.e. shallower) 
temperature gradients in the lower TR than those used by 
Shoub (1982,1983). This points to a potential problem with 
the formulation used here, which assumes a power law tail 
to exist at all temperatures. 

The low values of £ ta jj required to represent Shoub's 
(1982,1983) results are directly related to the steep temper- 
ature gradient in the low TR of the model atmosphere he 
used. Comparison of the atmospheric models used in the ra- 
diative transfer calculations presented here (VAL C, models 
S, X) with those used by Shoub (1982,1983) and Ljepojevic 
& Burgess (1990) suggest that values of £tail > 4.0 are more 
appropriate to realistic models. Given this argument, a pa- 
rameter space of 3.5 < £ ta jj < 5.0, 2.0 < r ta jj < 4.5 was 
explored in these calculations, with results for £tail > 4.0 
given greater credence. 

4.3 Collisional excitation and ionization rates 

In order to treat excitation and ionization by the non- 
Maxwellian EVDFs described above in the radiative trans- 
fer calculations, options were written into the multi code to 
compute the relevant rates. The collisional excitation and 
ionization rate coefficients for helium were calculated by in- 
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Figure 5. Dimcnsionlcss anglc-avcragcd EVDFs used to calculate 
enhanced collisional excitation rates in helium. (2-7rA;T/m) 3 / 2 /(£) 
is given as a function of £. Plotted in this way, the shape of 
the EVDF is independent of T. Dashed lines show Maxwellian 
EVDFs. The ranges of £tail anc ^ r tail use d i n the power law tail 
arc illustrated. 



tegrating electron collision cross-sections over the speed and 
the electron speed distribution. If the collisional excitation 
rate from level i to level j is Rij = NeC\j , the rate coefficient 
dj is given by 

poo 

dj = / 4nv 3 f(v)a(v)dv (3) 

where f(v) is the angle- averaged EVDF, a(v) is the speed- 
dependent collision cross-section for the transition, and Eij 
is the threshold energy of the transition. 

For an EVDF consisting of a Maxwellian distribution 
with temperature T for velocities smaller than v^, with a 
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power law decline of v p at higher velocities, dj becomes 




where the identities v-^ = £tail u th ana - P ~ (^/ r tail + ^) 
have been used to simplify the expression. If, at temperature 
T, vi < v(Eij) (i.e. the power law tail begins at an energy 
below threshold), then the first integral vanishes and the 
lower bound on the second integral becomes v(Eij). 

Analytical expressions for the collision cross-sections 
were used as the cross-sections were needed up to large elec- 
tron energies, where few numerical data exist for many of the 
transitions. Keeping the expressions for the rate coefficients 
in an analytical form also made it simple to make £tail ano - 
r tail f ree parameters, and kept the expressions transparent 
in their dependence on these parameters and on the tem- 
perature. The drawback of this approach is that for some 
transitions newer and better numerical data exist which are 
used in the standard atomic model described in Section 2, 
but not in these calculations of the rates due to non-local 
electrons. The approach taken by Anderson et al. (1996), in 
which the EVDF at each height is written as a weighted sum 
of Maxwellians, allowing use of numerical data tabulated for 
Maxwellian distributions, is an improvement in this respect, 
but the scarcity of high temperature data (for He I in par- 
ticular) is a problem for that method. 

The specific expressions for the rate coefficients depend 
on the cross-section used. For collisional ionization of He I 
and He II, expressions for the cross-sections given by Miha- 
las & Stone (1968) were used. Mihalas & Stone (1968) also 
provide bound-bound collision cross-sections for He II, using 
an easily-integrated semi-empirical formula due to Hinnov 
(1966), and for the optically allowed transitions in He I. For 
the forbidden transitions in He I, the expression due to Green 
(1966) given by Benson & Kulander (1972) was used. Benson 
& Kulander (1972) made order of magnitude estimates for / 
values used in the expression. I have made new estimates by 
matching as closely as possible the collision strengths given 
for T e = 6 x 10 3 K and T e = 2.5 x 10 4 K to more recently cal- 
culated numerical data used in the standard helium model 
(Lanzafame et al. 1993; Sawey & Berrington 1993). In cases 
where the integration above could not be performed an- 
alytically, it was performed numerically; such integrals are 
estimated to be correct to about 5 per cent. 

The modifications use collision rate data which differ 
from those in original helium model described in Section 2. 
When the modified rates are used with a very high value 
of £ ta y to approximate the Maxwellian regime, this results 
in the computed intensities of the He I resonance series and 
the He n 303. 8-A line differing from those calculated using 
the original model by at most 10 per cent. Owing to ap- 
proximations made for the forbidden transitions of He I the 
non-Maxwellian modifications cannot be used to investigate 
possible non-Maxwellian effects on the triplet lines of He I. 

Enhanced collision rates were not adopted for all tran- 
sitions in He I, as the suprathermal tail electrons will have 
little effect on transitions for which the threshold energy is 
much lower than the energy at which (locally) the power law 



tail of the EVDF begins. In such a case, even if the collision 
cross-section is larger at energies in the power law tail, the 
rapid decline of the EVDF above threshold (see Figure 5) 
means that excitation will be dominated by the Maxwellian 
part of the distribution. Hence the criterion for a transition 
to be relatively unaffected by the presence of an enhanced 
suprathermal tail of the type studied here is for its energy 
Eij to satisfy the inequality 

< al ?tail fcT ( 5 ) 

where (^g^kT ( = mv t/^ ^ s tlle ener gy at which the 
suprathermal tail departs from the Maxwellian distribution 
at temperature T. Even in the most extreme cases encoun- 
tered in the models tested here (i.e. T ~ 4000 K, £tail = 3.5), 
this inequality is satisfied for transitions in He I with An = 
except for those with n = 2, and for transitions between He I 
levels with n — 4 and n = 5. Collisional excitation rates for 
transitions satisfying the inequality are calculated assuming 
a Maxwellian distribution. 



5 RESULTS WITH NON-MAXWELLIAN 
COLLISION RATES 

The integrated intensities computed for the 303. 8-A, 584.3- 
A, and 537.0-A lines are given in Tables 3-6 for the VAL 
C model atmosphere and the network models S and X, for 
different values of the parameters £tail an< ^ r tail" P r °fil es °f 
the 584. 3-A and 303. 8-A lines computed for different values 
of £tail' w ith r tail ~ ^-5, are presented in Figures 6 and 7. 

Results for Ctail = ^-5 are given only where f+ a jj = 3.5 
and 3.0, in order to test the EVDFs computed by Shoub 
(1982,1983). This value of £tail ' s probably unrealistic in 
the lower TR, given the calculations of Ljepojevic & Burgess 
(1990), whose results are approximated here by r^ a jj = 2.0 
and 4.5, £tail — 4-0. The results for £tail ~ ^-5 should be 
regarded with caution, as the form of the modifications to 
the radiative transfer code results in there being significant 
effects in chromospheric regions of the model atmospheres, 
where in reality high energy electrons from the upper TR 
are unlikely to penetrate. Non-Maxwellian tails might exist 
at these heights if high energy electrons were accelerated 
locally by MHD processes, but without detailed models it is 
not known if the resultant distributions would resemble the 
ones used here. 

Even if the results for = 3.5 are discounted, com- 
parison of Tables 3-6 with Table 2 shows that collisional 
excitation by non-Maxwellian electron distributions can lead 
to significant enhancements of the helium resonance line in- 
tensities over those when Maxwellian distributions are as- 
sumed. The results also show that non-Maxwellian EVDFs 
could be present in the TR without producing a strong sig- 
nal in the helium lines (compare Tables 6 and 2). Calcula- 
tions for Ctail ~ 5-0, ?"tail = parameters characteristic of 
EVDFs computed by Ljepojevic & Burgess (1990), suggest 
that the only observable effect of such a distribution would 
be an increase in the intensity of the He I 537.0-A line with 
respect to the case of Maxwellian excitation. 

The trends in the calculated intensities are largely un- 
derstandable. For a given r^jj, a reduction of £tail generally 
causes increases in line intensities, as the power law tail be- 
comes significant at lower electron energies, and at a given 
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Table 3. Intensities and ratios of helium lines calculated using 
non-Maxwellian EVDFs with r tai [ = 4.5. 



Table 5. Intensities and ratios of helium lines calculated using 



non-Maxwellian EVDFs with r. 



tail 



3.0. 





?tail : 


5.0 


4.5 


4.0 




5tail : 


5.0 


4.5 


4.0 


3.5 


Model 


A 


Computed intensity 


Model 


A 




Computed intensity 


atmosphere 


(A) 


(erg . 


— 9 — 1 

cm s 


sr- 1 ) 


atmosphere 


(A) 


(erg cm 2 


s 1 sr 1 


) 


VAL C 


584.3 


400 


329 


675 


VAL C 


584.3 


411 


357 


639 


2243 




537.0 


51.9 


99.7 


192 




537.0 


44.7 


90.6 


179 


800 




303.8 


2026 


2032 


3259 




303.8 


2025 


2028 


2145 


12200 


Network S 


584.3 


148 


130 


349 


Network S 


584.3 


152 


139 


308 


1194 




537.0 


26.7 


48.4 


134 




537.0 


23.7 


43.8 


113 


561 




303.8 


1165 


1165 


1776 




303.8 


1165 


1165 


1227 


5709 


Network X 


584.3 


84.3 


85.3 


392 


Network X 


584.3 


86.1 


84.7 


313 


1527 




537.0 


11.4 


23.4 


129 




537.0 


9.87 


19.9 


93.7 


606 




303.8 


994 


999 


1792 




303.8 


994 


994 


1048 


6691 






/(A)//(584.3 A) 








/(A)//(584.3 A) 




VAL C 


537.0 


0.130 


0.303 


0.284 


VAL C 


537.0 


0.109 


0.254 


0.280 


0.357 




303.8 


5.07 


6.18 


4.83 




303.8 


4.93 


5.68 


3.36 


5.44 


Network S 


537.0 


0.180 


0.372 


0.384 


Network S 


537.0 


0.156 


0.315 


0.367 


0.470 




303.8 


7.87 


8.96 


5.09 




303.8 


7.66 


8.38 


3.98 


4.78 


Network X 


537.0 


0.135 


0.274 


0.329 


Network X 


537.0 


0.115 


0.235 


0.299 


0.397 




303.8 


11.8 


11.7 


4.57 




303.8 


11.5 


11.7 


3.35 


4.38 



Table 4. Intensities and ratios of helium lines calculated using 
non-Maxwellian EVDFs with r tai [ = 3.5. 



Table 6. Intensities and ratios of helium lines calculated using 
non-Maxwellian EVDFs with r tai [ = 2.0. 
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2025 
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584.3 
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Network S 


584.3 
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248 
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24.6 
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118 


479 




537.0 


22.0 


36.8 


115 
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584.3 


84.7 


83.6 


336 
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89.0 
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VAL C 


537.0 


0.113 


0.273 


0.266 


0.363 


VAL C 


537.0 


0.099 


0.200 


0.392 




303.8 


4.84 


5.81 


3.34 


15.2 




303.8 


4.94 


5.33 


3.56 


Network S 


537.0 


0.164 


0.339 


0.365 


0.437 


Network S 


537.0 


0.144 


0.250 


0.463 




303.8 


7.77 


8.63 


3.91 


12.5 




303.8 


7.61 


7.93 


4.75 


Network X 


537.0 


0.122 


0.251 


0.304 


0.406 


Network X 


537.0 


0.102 


0.187 


0.374 




303.8 


11.7 


11.9 


3.29 


13.1 




303.8 


11.2 


11.1 


4.28 



temperature increases the proportion of electrons with ener- 
gies high enough to collisionally excite a given line. A steeper 
power law (smaller r^ a ^) tends to suppress He II with respect 
to He I (compare Tables 3 and 6), as higher energy electrons 
are present in relatively smaller numbers. 

For any particular r ta jj, the He II resonance line inten- 



sity shows a sharper increase at low ^ a jj than do the He I 
lines, which is due to a significant change in the ionization 
balance. Enhanced line emission can occur, not just because 
of the increased collisional excitation of the relevant transi- 
tion, but also because the shape of the peak in the ionization 
fraction can be changed by the enhanced ionization rates at 
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low temperatures. Thus emission occurs over a larger, more 
dense region in He II, whereas the region of peak He I ioniza- 
tion fraction can only be narrowed. At the lower limit of the 
£tail tested, the ionization fractions at a given temperature 
depart from the standard model results by factors approach- 
ing, and in some places exceeding, an order of magnitude. 
As stated above, this may be unrealistic in the regions of 
He I line formation. 

The formation of the He n resonance line in higher den- 
sity regions when Ctail ' s sman can lead to its profile becom- 
ing self-reversed. For £tail > 4.0, the profile is approximately 
Gaussian, and very similar to that predicted by calculations 
made with the non-Maxwellian modifications. The compar- 
sion is shown for model S in Figure 7; results from the other 
models are similar. Figure 6 shows that the He I resonance 
line is self-reversed in all cases of non-Maxwellian excitation. 

The somewhat surprising reduction of the 584. 3-A 
line intensity (compared to calculations with 'normal' 
Maxwellian rates) for higher ^+ a j] (compare Tables 4 and 
2) appears to be due to a similar effect. The shift in the 
ionization balance suppresses the Maxwellian contribution 
to collisional excitation. The effect is smallest in model X, 
where Maxwellian collisional excitation is least significant 
because of the smaller amount of material in the relevant 
temperature range. The 584. 3-A line intensity can be re- 
duced if the tail has a value of £ tau too high for increased 
excitation by suprathermal electrons to compensate. Other 
lines of the He I resonance series do not show this, as the 
increased excitation has a greater effect than the shift in ion- 
ization (the suppressed Maxwellian contribution is relatively 
less significant). 

The effect of enhanced ionization on the ioniza- 
tion balance also complicates the trend in the ratio 
1(537.0 A)/7(584.3 A), which is expected to decrease with 
decreasing r ta y. For low values of £tail> tne trend can re- 
verse, as the change in the ionization balance affects the 
584.3-A line more than the 537.0-A line, allowing /(537.0 A) 
to grow more quickly with decreasing Ctail- 

5.1 Comparison with observations 

The results of the calculations may be compared with the 
mean line intensities observed by MJ99, which are given in 
Table 1. It can be seen that no particular combination of 
£^ a ;i and r^jj leads to a good match to observations in all 
three lines for any of the model atmospheres tested. 

In models in which the calculated 303. 8-A line inten- 
sity is comparable to observed values, the He I line ratio 
1(537.0 A)//(584.3 A) is much larger than observed (inde- 
pendent of the CDS calibration assumed) , and the computed 
303. 8- A line profile is in general self-reversed. Such reversals 
are not observed, but would probably not be resolved by the 
instruments on board SOHO (Smith 2000). Models showing 
enhancements of the 303. 8-A line intensity comparing well 
with observations have other problems. The values of Ctail 
required may be unrealistically low, but if such such a tail 
were present, significant collisional enhancement of the He n 
1640. 4-A multiplet is predicted. Intensities much larger than 
those observed, in some cases by orders of magnitude, are 
produced when £ ta y < 4.5 (independent of r ta y). This is in 
part due to the increased fraction of He n at lower temper- 
atures where recombination dominates the formation of the 
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Figure 6. He I 584.3-A line profiles produced by models VAL 
C (top), and network S (middle) and X (bottom), for non- 
Maxwellian EVDFs with r tail = 3.5 and £ tail = 5.0 (solid), 
?tail = 4 - 5 (dotted), § tail = 4.0 (dashed), £ tail = 3.5 (dot- 
dashed) . 

multiplet. Small departures from Maxwellian EVDFs might 
explain the Balmer a intensity, but results found here imply 
that non-Maxwellian collisional excitation could not be re- 
sponsible for the observed resonance line intensity without 
producing too much Balmer a emission. It remains to be 
seen whether non-Maxwellian collisional excitation is consis- 
tent with observations of the profiles and relative intensities 
of the components of the 1640. 4-A multiplet. An improved 
treatment of the energy levels of He II would allow an inves- 
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Figure 7. He II 303. 8-A line profiles produced using network 
model S, with Maxwcllian collision rates (top), and with a non- 
Maxwcllian EVDF with r tail = 3.5, Ctail = 3 - 5 (bottom). 



tigation of this point. The excitation of the 1640. 4- A multi- 
plet by non-local electrons will also be tested against stellar 
observations in future work. 

Trends in the results given in Tables 3-6 may be com- 
pared with observed variations in line intensities and ratios 
over different regions of the Sun. The effectiveness of the 
excitation mechanism depends on the temperature gradient 
and electron density, which have different values in the net- 
work and internetwork, in the quiet Sun and coronal holes. 
Given the significant disagreements between the calculations 
and observations described above, an approximate approach 
was followed in preference to constructing models of different 
parts of the atmosphere. The effects of different atmospheric 
parameters on the suprathermal electron tail are considered, 
and the results given in Tables 3-6 are used to infer how 
the helium line intensities might vary as a result. 

Changes in the suprathermal tail may be parametrized 
in simple terms as variations of Ctail ano - r tail' ^ ne var i~ 
ation of r^- a ji appears to depend on the energy balance in 
the atmosphere (Shoub 1982), and given the current lack 
of detailed understanding of the heating of the atmosphere, 
it is difficult to predict how conditions in different regions 
might affect r tai j. On the other hand, Shoub's (1982) BGK 
calculations provide an expression for the variation of £tail 
with easily understood physical quantities: 



where L is the thickness of the transition region and T n is 
the temperature at its upper edge; T^/L corresponds to the 
average temperature gradient. £ ta y is therefore decreased 
(and the influence of the suprathermal tail increased) by an 
increase in T n or the temperature gradient AT /Ah, or by a 
decrease in the electron pressure. This relation emerges from 
analytical calculations in a simplified system, but it shows 
expected dependences. The mean free path of an electron 
with a given velocity should increase with decreasing elec- 
tron density, and the change in temperature over that path 
will increase as the temperature gradient steepens. Thus 
such changes increase the numbers of electrons from a region 
at Ti reaching a region of lower temperature Ti. Increasing 
Ti shifts the peak of the distribution function to higher en- 
ergies, increasing the number of electrons in the tail with a 
given velocity. This also leads to an increase in the number 
of electrons with mean free paths large enough to influence 
lower temperature regions. 

The dependence of £t a il on eacn of the parameters is 
weak, but Tables 3-6 show that the helium line intensi- 
ties can be sensitive to relatively small changes in Ctaih ^ 
changes in £tail are assumed to control the variation of he- 
lium line emission with position on the Sun, the predicted 
variations in intensities may be tested against observations. 

With respect to average conditions in the quiet Sun, in 
coronal holes the electron pressure is 2 - 3 times smaller, 
the mean temperature gradient up to an order of magnitude 
smaller, and the coronal temperature at least 30 per cent 
smaller (Munro & Withbroe 1972; Jordan et al. 2001). Ac- 
cording to the expression in equation (6), £tail would then 
be a factor of up to ~ 1.5 larger in coronal holes than in the 
quiet Sun. Such a change spans the values investigated here; 
Tables 3-6 show that over such a range, very large changes 
in the helium line intensities can occur. The helium reso- 
nance line intensities are observed to be factors of 1.5 - 2.0 
smaller in coronal holes than in the quiet Sun (Peter 1999; 
Jordan et al. 2001). The present calculations suggest that 
such a decrease in He I line intensities would be accompanied 
by a large decrease in the 7(537.0 A)/7(584.3 A) line ratio 
(perhaps by a factor of two or more). This would outweigh 
the small increase in the ratio expected with the reduction 
of coronal illumination. Observations (Jordan et al. 2001) 
suggest that the ratio is slightly larger in coronal holes than 
in the quiet Sun, contrary to the predictions of models dom- 
inated by non-Maxwellian collisional excitation. Changes in 
£tail °f a f ac t° r °f 1-5 also generally lead to significant vari- 
ations in the calculated ratio of the He I and He 11 resonance 
lines. In contrast, observations show no consistent variation 
in this ratio between coronal holes and the quiet Sun; the 
mean observed ratio is almost exactly the same in the two 
regions (Jordan et al. 2001). 

The observed trends in line ratios in the quiet Sun may 
also be compared with the predictions of the calculations. 
In the quiet Sun MJ99 observed only small changes in the 
ratio 7(537.0 A)/7(584.3 A) for changes in 7(584.3 A) of a 
factor of 5, the ratio decreasing with increasing 1(584.3 A). 
The ratio 7(303.8 A)/7(584.3 A) is also observed to vary 
by relatively small amounts, and while consistent trends are 
not seen clearly, small values of the ratio often coincide with 
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large values of 7(584.3 A). For large variations in 1(584.3 A), 
the calculations predict 1(537.0 A)/7(584.3 A) to increase 
with increasing 7(584.3 A) in all but a small part of the 
parameter space, particularly for the (probably more realis- 
tic) smaller values of r^g^. The variation of the ratio is also 
larger than in observations. The predicted trend appears in- 
evitable given the exciting mechanism, since the power law 
tail provides relatively greater numbers of electrons (with re- 
spect to a Maxwellian distribution) at the excitation energy 
of the 537.0-A line than at the energy of the 584. 3-A line. 
The models also predict much larger variations in the ratio 
7(303.8 A)/7(584.3 A) than are observed for given changes 
in 7(584.3 A). Thus the predicted variations of the helium 
line ratios are not consistent with observations. 

The comparisons made so far with quiet Sun observa- 
tions have ignored the question of whether would be 
expected to vary in a manner which would produce the ob- 
served variations in intensities between cell interior and net- 
work regions. In considering the dependence of Ctail on ^ ne 
conditions in the quiet Sun transition region, the likely mag- 
netic field geometry should be taken into account together 
with the changes in temperature gradient and electron pres- 
sure. All calculations described above assume any magnetic 
field to be aligned with a radial temperature gradient. In 
general, electrons streaming from the corona and upper TR 
would tend to follow the magnetic field, which in the lower 
TR and chromosphere is concentrated in the network. The 
model of Gabriel (1976) suggests that the temperature gra- 
dient in the direction of the magnetic field is largest in the 
centre of the network and decreases towards the edges. The 
model also predicted lower pressures in the network, while 
MJ99 observed such a trend. Both factors would cause £tail 
to be smallest at the centres of network boundaries, sug- 
gesting that enhancement of the helium line intensities by 
non-local electrons would be concentrated in the centre of 
the network. Such behaviour is not observed; indeed the he- 
lium network appears to be unusually broad (e.g. Brueckner 
& Bartoe 1974; Gallagher et al. 1998), and the enhancement 
of the helium lines with respect to other transition region 
lines appears to increase towards cell interiors (MJ99). This 
observed behaviour could be due to radiative transfer ef- 
fects not considered in this simple argument, but it presents 
further problems for the proposed enhancement process. In 
contrast, it is argued in Paper I that enhancement by up- 
ward transport of helium ions in an expanding network field 
would lead naturally to the network appearing wider in the 
helium resonance lines than in other low TR lines. 

5.2 Approximations in the EVDF 

Given that the main approximation in the present work is 
that the form of the EVDF is assumed rather than calculated 
self-consistently, the extent to which the form chosen for 
the EVDF affects the computed line intensities should be 
considered. Two features of the EVDF in particular were 
investigated: the behaviour of the tail at very high velocities, 
and the dependence (or otherwise) of on temperature. 

Shoub's (1982) analytical results show a down-turn in 
the high energy tail, where the power law ceases to be a good 
description, at £ ~ 4.2 x 10 3 T -1 ^ 2 . The numerical calcula- 
tions do not show this feature, these values of £ being outside 
the range of calculation, but it is plausible that such a down- 



turn should exist, as the above limit on £ decreases with in- 
creasing T, so that the EVDF tends towards the Maxwellian 
distribution assumed to exist in the near-isothermal corona. 
When calculations were made with a cut-off in the EVDF 
at the above £(T), significant effects on the helium reso- 
nance lines were found for small Ctaik ^ ne decreased extent 
of the power law tail at higher temperatures reduces ioniza- 
tion rates from He I to He II, decreasing the extent of the 
region where He II dominates the ionization balance. Con- 
sequently, the intensity of the 303. 8-A line is decreased by 
up to a factor of 2, and the intensities of the He I lines are 
increased by up to 50 per cent. The cut-off is a crude device 
and calculations with ?^ a ;j < 3.5 may well approximate the 
effects of the down-turn more accurately, but the sensitiv- 
ity revealed here suggests that this part of the distribution 
may need closer attention in the future. The present results 
imply that including the down-turn properly in calculations 
would not improve the fit to observations. 

The use in this study of EVDFs with Ctail independent 
of T is an approximation suggested by the general form of 
solutions to the Boltzmann equation. Shoub's (1982) analyt- 
ical solutions suggest £tail ~ J 11 / 12 (thus it changes by only 
20 per cent for an order of magnitude change in T) . Includ- 
ing such a dependence in calculations would be expected to 
worsen agreement with observations. The observed He II res- 
onance line intensity can only be produced by low values of 
£taib kut t ne ratio of the 584. 3-A and 537.0-A lines favours 
higher values, which would require Ctail ^° var y with T in 
the opposite sense. Ljepojevic & Burgess' (1990) results do 
show some increase in for the lowest temperatures for 
which they present the EVDF. If were lower at the 

higher T of He II line formation than at the lower T of He I 
line formation, a better fit might be obtained, although the 
problem with the He II Balmer a intensity would remain. 
This problem is certainly seen in the results of calculations 
made with enhanced collision rates operating only in He II 
and not in He I, approximating very roughly the situation 
in which is increased at low temperatures. 



6 SUMMARY AND CONCLUSIONS 

A detailed study has been made of the formation of the 
helium resonance lines, involving radiative transfer calcula- 
tions with a 36 level model helium atom. Calculations using 
the the VAL C (average quiet Sun) model atmosphere, with 
a coronal radiation field characteristic of the quiet Sun based 
on the solar irradiance model of Tobiska (1991), come close 
to reproducing the observed He I 584. 3-A line intensity, but 
produce intensities in the He I 537.0-A and He II 303. 8-A 
lines smaller than are observed. The predicted self-reversal 
of the 584. 3-A line is not observed, but it is shown that it 
may be largely obscured by instrumental effects. Calcula- 
tions using the VAL D (network) model produce results in 
better agreement with observed network intensities, but the 
VAL models both produce poor matches to other transition 
region lines (e.g. of C II and Si in). 

New model atmospheres (S and X) appropriate to net- 
work boundary regions were constructed to be more con- 
sistent with observed intensities of TR lines of species other 
than helium. Calculations using these new models with quiet 
coronal illumination fail to reproduce the observed quiet Sun 
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network He n resonance line intensity by factors of 7-8 (3- 

4 when the most recent CDS calibration is used). Model X, 
based on the £ Boo A EMD, produces an intensity for the 
He I resonance line a factor of 7.5 smaller than observed, 
but model S, based on the MJ99 solar network EMD, gives 
better agreement. The computed intensity for the 584.3 A 
line is a factor of 4 too small, but the computed line profile is 
closest to the observed one, having only a small self-reversal 
that could be completely obscured by instrumental broaden- 
ing. The ratio 7(537.0 A)/7(584.3 A) computed using model 

5 also matches observations most closely. 

The optical thickness of the He I and He n resonance 
lines may explain the unusual width of the network observed 
in those lines if photons are scattered from the edges of the 
network. Studies of two-component model atmospheres are 
required to investigate this point. 

In all of the model atmospheres the formation of the 
He II resonance line is dominated by collisional excitation, 
but recombination is more important in the formation of the 
He II 1640. 4-A multiplet. Recombination also contributes 
significantly to the intensity of the He I 584. 3-A line, but 
does not dominate its formation except in model X. Direct 
collisional excitation from the ground and indirect excita- 
tion by radiative transitions from other collisionally excited 
singlet levels also contribute significantly to the formation 
of the line. Similar formation processes are important in the 
He I 537.0-A line. 

Increasing the intensity of the coronal radiation field 
causes increases in the computed He I line intensities and 
generally produces a decrease in the 7(537.0 A)/7(584.3 A) 
ratio. The He II resonance line intensity is generally reduced 
by increased coronal illumination, as a shift in the peak 
He II ionization fraction to lower temperatures reduces col- 
lisional excitation rates. Removing the coronal illumination 
produces very little change in the computed He II resonance 
line intensity. An absence of coronal radiation does cause a 
reduction in the computed intensities of the He I lines, an 
effect approaching the observed reduction in model S, ex- 
ceeding it in model X. The line ratio 7(537.0 A)/7(584.3 A) 
is observed to be slightly larger in coronal holes than in 
the quiet Sun, and this feature is reproduced by the models 
with zero coronal illumination. The quite different responses 
to changes in coronal irradiance of the He I and He II res- 
onance lines suggest that this is not the factor controlling 
factor the intensities of both lines, given that they show sim- 
ilar behaviour, both in the quiet Sun and coronal holes, and 
that the ratio of their observed intensities does not change 
by large amounts (MJ99, Jordan et al. 2001). 

Given that the intensities of the helium lines with re- 
spect to other TR lines are not explained by the 'standard' 
models used in the present work, modifications were made 
to the radiative transfer code to allow the effects of non- 
Maxwellian EVDFs on collision rates to be included in calcu- 
lations of the helium line intensities. The EVDFs tested are 
of a form postulated to exist in the transition region because 
of the extreme temperature gradient (e.g. Shoub 1983); they 
are Maxwellian at thermal velocities, with a power law tail 
at suprathermal velocities. The parameters of the power law 
tail were varied in a range suggested by the results of Shoub 
(1982,1983) and Ljepojevic & Burgess (1990), allowing an 
exploration of EVDFs that might plausibly be present in the 
model atmospheres studied here. 



It is found that this simulation of the presence of an in- 
creased population of suprathermal electrons in the lower 
transition region can lead to significant enhancements of 
the helium resonance line intensities compared with mod- 
els assuming Maxwellian electron distributions. None of the 
models tested, however, simultaneously reproduce the ob- 
served intensities of the three lines studied. Calculations in 
which the He II resonance line intensity is approximately 
reproduced also produce a much larger intensity in the 
He II 1640. 4-A multiplet than is observed. The enhanced 
collision rates produce large values of the He I line ratio 
7(537.0 A)/7(584.3 A), as the increase over Maxwellian rates 
is relatively greater in the higher energy transition. In some 
cases, where the power law tail is relatively insignificant, 
this produces better matches between the computed and 
observed ratios, but in cases where the suprathermal tail 
dominates excitation the computed ratio is generally much 
larger than observed. The calculations also predict that the 
ratio would increase with the absolute intensities of the lines, 
which is the opposite of the observed trend (MJ99). 

A consideration of how local atmospheric conditions 
might affect the EVDF at a given temperature in the low 
TR suggests that a reduction of the helium resonance line 
intensities in coronal holes with respect to the quiet Sun 
would be expected if excitation by non-local electrons dom- 
inates line formation. However, the models also predict that 
a decrease in absolute He I intensities of the observed or- 
der would be accompanied by a decrease in the line ratio 
7(537.0 A)/7(584.3 A); the opposite is observed (Jordan 
et al. 2001). This excitation process would also probably pro- 
duce a larger variation in the ratio 7(303.8 A)/7(584.3 A) 
than is observed between coronal holes and the quiet Sun 
(Jordan et al. 2001). Similar considerations suggest that the 
appearance of the network in the helium lines is not easily 
explained in the case of excitation by non-local electrons. 

In general, the calculations predict that significant de- 
partures from Maxwellian EVDFs of the form investigated 
would produce signatures in the helium line ratios which 
contradict observations. Small departures from Maxwellian 
EVDFs are not ruled out, but I conclude that collisional 
excitation by suprathermal electrons in EVDFs of the form 
examined here is unlikely to dominate the formation of the 
helium resonance lines. Enhancement by non-thermal trans- 
port of helium atoms and ions, as investigated in Paper I, 
appears to be a more promising explanation of the helium 
resonance line intensities. 

The calculations performed in the present work contain 
a number of approximations, and the above conclusions are 
therefore contingent on a more sophisticated study of the 
excitation process. Ideally, such an investigation would be 
totally self-consistent, solving the Boltzmann equation for 
the EVDF and the equation of radiative transfer simulta- 
neously. Models in which suprathermal electron populations 
are generated by reconnection events or MHD processes in 
the low TR (rather than by the steep TR temperature gra- 
dient) also deserve examination. The apparent sensitivity of 
the He II 1640.4 A multiplet to the presence of suprathermal 
electrons (in part through the altered ionization balance) 
suggests that these lines also deserve further investigation. 
In particular, the plausibility of enhanced suprathermal ex- 
citation could be tested against the observed profile of the 
multiplet and also against stellar observations. 
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APPENDIX 

The new model atmospheres described in Section 2.2 arc 
presented in the following pages. The tables give the loga- 
rithm (to base 10) of the mass column density m (in units 
of g cm -2 ), the electron temperature, the neutral hydrogen, 
proton, and electron number densities, the microturbulent 
velocity « T , and the height, h, above optical depth unity at 
5000 A. 
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Table 7. Network model S 



log m T e (K) N[j (cm -3 ) N p (cm -3 ) N e (cm -3 ) v T (km s -1 ) h (km) 
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log m T e (K) N H (cm- 3 ) N p (cm- 3 ) iV e (cm" 3 ) v T (km s" 1 ) h (km) 



-5 


11289 


14231.2 


1 


021 


+10) 


3 


428 


+10) 


3 


872 


'+10) 


8.741 


2 


04160 


;+03) 


-5 


11272 


13448.8 


1 


201 


+10) 


3 


508 


+10) 


3 


973 


'+10) 


8.736 


2 


04133 


;+03) 


-5 


11269 


13316.1 


1 


234 


+10) 


3 


523 


+10) 


3 


994 


;+io) 


8.735 


2 


04128 


;+03) 


-5 


11260 


12860.5 


1 


348 


+10) 


3 


574 


+10) 


4 


059 


;+io) 


8.733 


2 


04113 


;+03) 


-5 


11250 


12467.5 


1 


454 


+10) 


3 


621 


+10) 


4 


113 


'+10) 


8.731 


2 


04098 


;+03) 


-5 


11240 


12012.6 


1 


582 


+10) 


3 


682 


+10) 


4 


174 


;+io) 


8.728 


2 


04084 


;+03) 


-5 


11230 


11555.1 


1 


718 


+10) 


3 


759 


+10) 


4 


220 


'+10) 


8.725 


2 


04070 


;+03) 


-5 


11220 


11097.5 


1 


861 


+10) 


3 


862 


+10) 


4 


233 


;+io) 


8.723 


2 


04057 


'+03) 


-5 


11210 


10813.9 


1 


958 


+10) 


3 


937 


+10) 


4 


220 


'+10) 


8.720 


2 


04044 


;+03) 


-5 


11191 


10700.7 


2 


016 


+10) 


3 


959 


+10) 


4 


203 


;+io) 


8.719 


2 


04019 


;+03) 


-5 


11171 


10599.0 


2 


071 


+10) 


3 


978 


+10) 


4 


188 


'+10) 


8.718 


2 


03995 


;+03) 


-5 


11132 


10405.0 


2 


176 


+10) 


4 


012 


+10) 


4 


164 


;+io) 


8.717 


2 


03947 


;+03) 


-5 


11053 


10007.5 


2 


393 


+10) 


4 


075 


+10) 


4 


142 


;+io) 


8.713 


2 


03853 


;+03) 


-5 


10975 


9610.0 


2 


620 


+10) 


4 


123 


+10) 


4 


151 


;+io) 


8.710 


2 


03763 


;+03) 


-5 


10881 


9529.3 


2 


727 


+10) 


4 


107 


+10) 


4 


131 


;+io) 


8.702 


2 


03659 


'+03) 


-5 


10788 


9450.1 


2 


830 


+10) 


4 


093 


+10) 


4 


114 


;+io) 


8.694 


2 


03556 


'+03) 


-5 


10600 


9288.7 


3 


029 


+10) 


4 


072 


+10) 


4 


088 


'+10) 


8.678 


2 


03352 


;+03) 


-5 


10413 


9128.8 


3 


225 


+10) 


4 


054 


+10) 


4 


068 


'+10) 


8.661 


2 


03153 


;+03) 


-5 


10126 


8883.4 


3 


522 


+10) 


4 


034 


+10) 


4 


046 


;+io) 


8.637 


2 


02856 


;+03) 


-5 


09864 


8658.3 


3 


793 


+10) 


4 


022 


+10) 


4 


032 


;+io) 


8.614 


2 


02592 


;+03) 


-5 


09502 


8480.8 


4 


077 


+10) 


3 


996 


+10) 


4 


006 


;+io) 


8.593 


2 


02237 


;+03) 


-5 


08913 


8365.4 


4 


402 


+10) 


3 


949 


+10) 


3 


960 


;+io) 


8.571 


2 


01672 


'+03) 


-5 


08323 


8250.0 


4 


721 


+10) 


3 


912 


+10) 


3 


922 


'+10) 


8.550 


2 


01118 


;+03) 


-5 


07604 


8125.0 


5 


094 


+10) 


3 


872 


+10) 


3 


883 


'+10) 


8.525 


2 


00456 


;+03) 


-5 


06885 


8000.0 


5 


468 


+10) 


3 


839 


+10) 


3 


850 


;+io) 


8.500 


1 


99807 


'+03) 


-5 


05586 


7870.0 


6 


057 


+10) 


3 


787 


+10) 


3 


799 


;+io) 


8.460 


1 


98663 


;+03) 


-5 


04287 


7740.0 


6 


649 


+10) 


3 


749 


+10) 


3 


762 


;+io) 


8.420 


1 


97548 


'+03) 


-5 


02682 


7635.0 


7 


371 


+10) 


3 


729 


+10) 


3 


742 


;+io) 


8.320 


1 


96206 


;+03) 


-5 


01077 


7530.0 


8 


108 


+10) 


3 


726 


+10) 


3 


739 


;+io) 


8.220 


1 


94901 


'+03) 


-4 


98549 


7420.3 


9 


204 


+10) 


3 


745 


+10) 


3 


760 


;+io) 


8.099 


1 


92906 


;+03) 


-4 


96121 


7322.3 


1 


031 


+11) 


3 


787 


+10) 


3 


804 


;+io) 


7.976 


1 


91048 


;+03) 


-4 


93694 


7247.5 


1 


147 


+11) 


3 


850 


+10) 


3 


867 


;+io) 


7.834 


1 


89242 


'+03) 


-4 


88839 


7116.9 


1 


399 


+11) 


4 


006 


+10) 


4 


026 


;+io) 


7.577 


1 


85771 


'+03) 


-4 


82271 


7005.8 


1 


773 


+11) 


4 


242 


+10) 


4 


265 


;+io) 


7.323 


1 


81306 


;+03) 


-4 


76027 


6900.2 


2 


207 


+11) 


4 


477 


+10) 


4 


503 


'+10) 


7.081 


1 


77253 


;+03) 


-4 


66481 


6786.4 


3 


030 


+11) 


4 


839 


+10) 


4 


871 


;+io) 


6.733 


1 


71400 


;+03) 


-4 


55600 


6698.4 


4 


275 


+11) 


5 


215 


+10) 


5 


260 


;+io) 


6.357 


1 


65133 


;+03) 


-4 


43506 


6600.7 


6 


211 


+11) 


5 


566 


+10) 


5 


618 


;+io) 


5.938 


1 


58601 


;+03) 


-4 


30139 


6535.7 


9 


185 


+11) 


5 


898 


+10) 


5 


970 


;+io) 


5.523 


1 


51842 


'+03) 


-4 


15567 


6475.9 


1 


385 


+ 12) 


6 


294 


+10) 


6 


392 


;+io) 


5.096 


1 


44897 


;+03) 


-4 


07753 


6443.9 


1 


718 


+ 12) 


6 


524 


+10) 


6 


640 


'+10) 


4.877 


1 


41292 


;+03) 


-3 


99939 


6411.8 


2 


125 


+12) 


6 


812 


+10) 


6 


950 


;+io) 


4.659 


1 


37785 


;+03) 


-3 


91534 


6380.7 


2 


661 


+12) 


7 


222 


+10) 


7 


387 


;+io) 


4.432 


1 


34114 


;+03) 


-3 


83129 


6355.5 


3 


320 


+12) 


7 


739 


+10) 


7 


937 


;+io) 


4.219 


1 


30536 


;+03) 


-3 


65298 


6256.8 


5 


324 


+12) 


8 


481 


+10) 


8 


769 


;+io) 


3.777 


1 


23347 


;+03) 


-3 


47354 


6127.5 


8 


547 


+12) 


8 


741 


+10) 


9 


159 


;+io) 


3.356 


1 


16498 


'+03) 


-3 


30000 


6135.4 


1 


313 


+13) 


1 


125 


+11) 


1 


187 




2.960 


1 


10117 


;+03) 


-3 


13428 


6112.0 


1 


982 


+13) 


1 


356 


+11) 


1 


445 




2.592 


1 


04186 


;+03) 


-2 


97768 


6044.7 


2 


936 


+13) 


1 


489 


+11) 


1 


612 


'+11) 


2.262 


9 


O TO o o 

87388 


'+02) 


-2 


83040 


5955.0 


4 


258 


+13) 


1 


539 


+11) 


1 


702 


:+n) 


1.961 


9 


37643 


'+02) 


-2 


68992 


5851.9 


6 


071 


+13) 


1 


520 


+11) 


1 


729 


:+n) 


1.689 


8 


91579 


;+02) 


-2 


55407 


5736.4 


8 


529 


+13) 


1 


433 


+11) 


1 


690 


:+n) 


1.538 


8 


48232 


'+02) 


-2 


42406 


5589.2 


1 


187 


+14) 


1 


230 


+11) 


1 


527 


;+n) 


1.416 


8 


07900 


'+02) 


-2 


29839 


5444.7 


1 


635 


+14) 


1 


027 


+11) 


1 


365 


:+n) 


1.299 


7 


70060 


'+02) 


-2 


17917 


5291.1 


2 


222 


+14) 


8 


018 


+10) 


1 


175 


;+n) 


1.193 


7 


35239 


;+02) 


-2 


07437 


5132.7 


2 


923 


+ 14) 


5 


753 


+10) 


9 


920 


;+io) 


1.108 


7 


05562 


'+02) 


-2 


02936 


5061.5 


3 


291 


+ 14) 


4 


829 


+10) 


9 


219 


;+io) 


1.073 


6 


93045 


;+02) 
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log m 


T e (K) 


N H (cm- 3 ) 


N p (cm- 3 ) 


N e (cm" 3 ) 


t>T (km s ) 


h (km) 


-1.98435 


4987.5 


3.708 


+ 14) 


3.945 


+10) 


8.598 


+10) 


1.041 


6.80717(+02) 


-1.90597 


4858.6 


4.566 


+ 14) 


2.616 


+10) 


7.875 


+10) 


0.985 


6.59769(+02) 


-1.83566 


4746.6 


5.502 


+ 14) 


1.718 


+10) 


7.715 


+10) 


0.937 


6.41445(+02) 


-1.77067 


4646.6 


6.534 


+ 14) 


1.120 


+10) 


7.977 


+10) 


0.894 


6.24892(+02) 


-1.70766 


4549.5 


7.723 


+ 14) 


7.059 


+09) 


8.551 


+10) 


0.851 


6.09192(+02) 


-1.64196 


4469.2 


9.155 


+ 14) 


4.747 


+09) 


9.481 


+10) 


0.809 


5.93155(+02) 


-1.56959 


4402.8 


1.099 


+ 15) 


3.421 


+09) 


1.076 


+11) 


0.763 


5.75802(+02) 


-1.48889 


4328.7 


1.348 


+ 15) 


2.345 


+09) 


1.231 


+11) 


0.713 


5.56793(+02) 


-1.39929 


4297.7 


1.670 


+ 15) 


2.089 


+09) 


1.450 


+11) 


0.669 


5.35970(+02) 


-1.30098 


4280.3 


2.105 


+ 15) 


2.023 


+09) 


1.741 


+11) 


0.626 


5.13289(+02) 


-1.18674 


4277.6 


2.742 


+ 15) 


2.147 


+09) 


2.167 


+11) 


0.582 


4.87053(+02) 


-1.05256 


4295.8 


3.722 


+15) 


2.574 


+09) 


2.837 


+11) 


0.540 


4.56269(+02) 


-0.90681 


4354.9 


5.137 


+15) 


3.845 


+09) 


3.911 


+11) 


0.527 


4.22575(+02) 


-0.76297 


4423.9 


7.043 


+15) 


5.907 


+09) 


5.405 


+11) 


0.524 


3.88815(+02) 


-0.62093 


4492.0 


9.620 


+15) 


8.866 


+09) 


7.427 


+11) 


0.520 


3.54964(+02) 


-0.48037 


4586.4 


1.302 


+16) 


1.511 


+10) 


1.034 


+12) 


0.555 


3.20827(+02) 


-0.34126 


4683.4 


1.755 


+16) 


2.612 


+10) 


1.434 


+12) 


0.593 


2.86306(+02) 


-0.20325 


4782.3 


2.359 


+16) 


4.641 


+10) 


1.981 


+12) 


0.637 


2.51307(+02) 


-0.06709 


4930.9 


3.122 


+16) 


1.088 


+11) 


2.798 


+ 12) 


0.774 


2.15775(+02) 


0.03467 


5042.0 


3.852 


+16) 


2.104 


+11) 


3.626 


+ 12) 


0.877 


1.88360(+02) 


0.11629 


5131.1 


4.559 


+16) 


3.582 


+11) 


4.480 


+ 12) 


0.959 


1.65848(+02) 


0.18791 


5225.7 


5.268 


+16) 


6.114 


+11) 


5.477 


+ 12) 


1.037 


1.45685(+02) 


0.25064 


5327.2 


5.959 


+16) 


1.043 


+ 12) 


6.671 


+ 12) 


1.111 


1.27645(+02) 


0.30637 


5417.4 


6.650 


+16) 


1.647 


+ 12) 


8.034 


+ 12) 


1.176 


1.11293(+02) 


0.35814 


5527.4 


7.330 


+16) 


2.721 


+ 12) 


9.911 


+ 12) 


1.240 


9.57905(+01) 


0.41162 


5658.4 


8.083 


+16) 


4.692 


+ 12) 


1.281 


+ 13) 


1.307 


7.93916(+01) 


0.46173 


5781.2 


8.863 


+16) 


7.524 


+ 12) 


1.659 


+ 13) 


1.370 


6.36483(+01) 


0.50584 


5919.9 


9.566 


+16) 


1.208 


+ 13) 


2.209 


+13) 


1.428 


4.94471(+01) 


0.54316 


6068.0 


1.016 


+17) 


1.897 


+ 13) 


2.987 


+13) 


1.479 


3.71186(+01) 


0.57682 


6201.6 


1.073 


+17) 


2.766 


+13) 


3.944 


+13) 


1.525 


2.57236(+01) 


0.61048 


6335.2 


1.133 


+17) 


3.936 


+13) 


5.205 


+13) 


1.571 


1.40667(+01) 


0.64141 


6490.1 


1.186 


+17) 


5.715 


+13) 


7.088 


+13) 


1.614 


3.09627(+00) 


0.65687 


6603.4 


1.208 


+17) 


7.334 


+13) 


8.765 


+13) 


1.637 


-2.50818(+00) 


0.67234 


6716.7 


1.229 


+17) 


9.314 


+13) 


1.081 


+14) 


1.661 


-8.21229(+00) 


0.69582 


6888.8 


1.264 


+17) 


1.315 


+14) 


1.475 


+ 14) 


1.696 


-1.70655(+01) 


0.70757 


7011.2 


1.275 


+17) 


1.651 


+ 14) 


1.817 


+ 14) 


1.709 


-2.15903(+01) 


0.71931 


7145.4 


1.285 


+17) 


2.095 


+ 14) 


2.269 


+ 14) 


1.720 


-2.61991(+01) 


0.73755 


7353.9 


1.301 


+17) 


2.977 


+ 14) 


3.164 


+ 14) 


1.738 


-3.35360(+01) 


0.75580 


7562.4 


1.318 


+17) 


4.146 


+ 14) 


4.346 


+ 14) 


1.756 


-4.10870(+01) 


0.77133 


7754.9 


1.330 


+17) 


5.52E 


(+14) 


5.744 


+ 14) 


1.768 


-4.76907(+01) 


0.78686 


7953.0 


1.341 


+17) 


7.323 


+ 14) 


7.554 


+ 14) 


1.779 


-5.44673(+01) 


0.80107 


8134.1 


1.351 


+17) 


9.357 


+ 14) 


9.603 


+ 14) 


1.790 


-6.08198(+01) 


0.81517 


8313.9 


1.362 


+17) 


1.181 


+ 15) 


1.207 


+ 15) 


1.800 


-6.72742(+01) 
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Table 8. Network model X 



log m T e (K) N[j (cm -3 ) N p (cm -3 ) N e (cm -3 ) v T (km s -1 ) h (km) 



-5 


19116 


1584822.0 


5 


146 


+00) 


3 


474 


+08) 


4 


169 


+08) 


24.997 


1 


07156 


+04) 


-5 


17478 


1415053.0 


7 


568 


+00) 


4 


034 


+08) 


4 


840 


+08) 


24.997 


7 


92667 


+03) 


-5 


16217 


1261894.6 


1 


108 


+01) 


4 


648 


+08) 


5 


577 


+08) 


24.997 


6 


00867 


+03) 


-5 


15270 


1125551.6 


1 


613 


+01) 


5 


314 


+08) 


6 


376 


+08) 


24.997 


4 


72109 


+03) 


-5 


14560 


1004223.0 


2 


334 


+01) 


6 


039 


+08) 


7 


246 


+08) 


24.997 


3 


85669 


+03) 


-5 


14027 


896382.5 


3 


365 


+01) 


6 


831 


+08) 


8 


196 


+08) 


24.997 


3 


27644 


+03) 


-5 


13627 


800258.4 


4 


840 


+01) 


7 


699 


+08) 


9 


237 


+08) 


24.997 


2 


88681 


+03) 


-5 


13327 


715009.6 


6 


941 


+01) 


8 


647 


+08) 


1 


038 


+09) 


24.997 


2 


62507 


+03) 


-5 


13103 


639581.0 


9 


920 


+01) 


9 


681 


+08) 


1 


162 


+09) 


24.997 


2 


44905 


+03) 


-5 


12935 


572675.9 


1 


414 


+02) 


1 


081 


+09) 


1 


297 


+09) 


24.997 


2 


33070 


+03) 


-5 


12809 


513098.5 


2 


015 


+02) 


1 


204 


+09) 


1 


445 


+09) 


24.997 


2 


25091 


+03) 


-5 


12714 


461102.6 


2 


848 


+02) 


1 


337 


+09) 


1 


604 


+09) 


24.997 


2 


19713 


+03) 


-5 


12644 


415441.6 


4 


001 


+02) 


1 


478 


+09) 


1 


774 


+09) 


24.997 


2 


16072 


+03) 


-5 


12591 


375127.7 


5 


594 


+02) 


1 


630 


+09) 


1 


955 


+09) 


24.997 


2 


13602 


+03) 


-5 


12552 


341028.9 


7 


667 


+02) 


1 


784 


+09) 


2 


140 


+09) 


24.997 


2 


11918 


+03) 


-5 


12522 


310784.0 


1 


045 


+03) 


1 


947 


+09) 


2 


336 


+09) 


24.997 


2 


10766 


+03) 


-5 


12500 


285294.0 


1 


395 


+03) 


2 


109 


+09) 


2 


530 


+09) 


24.997 


2 


09975 


+03) 


-5 


12484 


263687.9 


1 


823 


+03) 


2 


268 


+09) 


2 


722 


+09) 


24.997 


2 


09426 


+03) 


-5 


12472 


244839.8 


2 


350 


+03) 


2 


428 


+09) 


2 


914 


+09) 


24.997 


2 


09044 


+03) 


-5 


12462 


229435.0 


2 


943 


+03) 


2 


577 


+09) 


3 


092 


+09) 


24.997 


2 


08776 


+03) 


-5 


12455 


214725.0 


3 


710 


+03) 


2 


737 


+09) 


3 


284 


+09) 


24.997 


2 


08561 


+03) 


-5 


12448 


201497.3 


4 


643 


+03) 


2 


898 


+09) 


3 


477 


+09) 


24.997 


2 


08378 


+03) 


-5 


12440 


185830.1 


6 


224 


+03) 


3 


129 


+09) 


3 


754 


+09) 


24.528 


2 


08206 


+03) 


-5 


12433 


168687.5 
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507 
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7 
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+09) 
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-5 


12386 
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3 


548 


+05) 


7 


885 


+09) 


9 


461 


+09) 
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07363 


+03) 
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12379 
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5 
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+05) 


8 


488 


+09) 


1 


018 
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2 


07298 
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-5 


12372 
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6 
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8 
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+09) 
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1 
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12354 
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525 


+06) 


1 


033 
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240 
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-5 
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+06) 
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110 
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1 


332 
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07048 
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-5 


12338 
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3 


648 


+06) 


1 


197 


+10) 


1 


437 
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2 


06997 


+03) 


-5 


12330 


41992.1 


6 


031 


+06) 


1 


297 


+10) 


1 


556 
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14.988 


2 


06949 


+03) 


-5 


12321 


38304.7 


1 
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+07) 


1 


421 


+10) 


1 


705 
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06904 
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-5 


12312 
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2 
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+07) 


1 


566 


+10) 


1 


880 
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2 


06860 


+03) 


-5 


12306 
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3 


630 


+07) 


1 


680 


+10) 


2 


017 


+ 10) 
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2 


06832 


+03) 


-5 


12300 


29930.5 


6 


410 


+07) 


1 


810 


+10) 


2 


172 


+ 10) 


12.805 


2 


06806 


+03) 


-5 


12294 


27636.5 


1 


207 


+08) 


1 


958 


+10) 


2 


348 


+ 10) 


12.310 


2 


06783 


+03) 


-5 


12288 


25957.6 


2 


029 


+08) 


2 


087 


+10) 


2 


491 


+ 10) 


11.891 


2 


06761 


+03) 


-5 


12282 


24080.4 


3 


888 


+08) 


2 


269 


+10) 


2 


635 


+ 10) 


11.423 


2 


06740 


+03) 


-5 


12275 
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6 


398 


+08) 


2 


422 


+10) 


2 


727 


+ 10) 


11.055 
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+03) 


-5 


12269 
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1 


025 


+09) 


2 
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+10) 


2 


835 
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+03) 


-5 


12260 
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+09) 


2 


720 


+10) 


3 


016 
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+03) 


-5 


12250 
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2 


759 


+09) 


2 


801 


+10) 


3 
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+ 10) 


9.954 


2 
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+03) 


-5 


12230 


17636.9 


4 
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+09) 


2 


947 


+10) 


3 


292 
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9.432 
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-5 


12211 


17553.1 


4 
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+09) 


2 


973 


+10) 


3 


319 
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9.413 


2 


06569 


+03) 


-5 
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4 


507 


+09) 


2 


994 


+10) 


3 


340 


+ 10) 


9.399 


2 


06527 


+03) 


-5 


12165 
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4 


428 


+09) 


3 
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+10) 


3 


366 


+ 10) 


9.380 


2 


06472 


+03) 


-5 


12139 
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4 


391 


+09) 


3 


039 


+10) 


3 


389 


+ 10) 


9.361 


2 


06417 


+03) 


-5 


12078 
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4 


398 


+09) 


3 


085 


+10) 


3 


439 


+ 10) 


9.316 


2 


06287 


+03) 


-5 


12013 
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4 


494 


+09) 


3 


129 


+10) 


3 


488 


+ 10) 


9.269 


2 
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+03) 


-5 
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4 
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+09) 


3 
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+10) 


3 


528 


+ 10) 


9.229 


2 
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-5 


11908 


16463.1 


4 


756 


+09) 


3 


195 


+10) 


3 


563 


'+10) 


9.193 


2 


05939 


;+03) 


-5 


11815 


16125.6 


5 


086 


+09) 


3 


251 


+10) 


3 


629 


'+10) 


9.124 


2 


05752 


;+03) 


-5 


11726 


15805.2 


5 


471 


+09) 


3 


304 


+10) 


3 


690 


;+io) 


9.059 


2 


05579 


;+03) 


-5 


11595 


15335.3 


6 


176 


+09) 


3 


377 


+10) 


3 


777 


;+io) 


8.964 


2 


05331 


;+03) 


-5 


11431 


14743.6 


7 


366 


+09) 


3 


460 


+10) 


3 


879 


;+io) 


8.844 


2 


05032 


;+03) 


-5 


11316 


14329.0 


8 


504 


+09) 


3 


506 


+10) 


3 


941 


;+io) 


8.761 


2 


04830 


;+03) 


-5 


11302 


14278.7 


8 


671 


+09) 


3 


511 


+10) 


3 


947 


;+io) 


8.750 


2 


04807 


;+03) 


-5 


11289 


14231.2 


8 


839 


+09) 


3 


514 


+10) 


3 


952 


;+io) 


8.741 


2 


04784 


'+03) 


-5 


11272 


13426.1 


1 


040 


+10) 


3 


615 


+10) 


4 


078 


'+10) 


8.736 


2 


04756 


;+03) 


-5 


11269 


13299.7 


1 


066 


+10) 


3 


632 


+10) 


4 


098 


;+io) 


8.735 


2 


04751 


;+03) 


-5 


11260 


12840.8 


1 


164 


+10) 


3 


695 


+10) 


4 


173 


'+10) 


8.733 


2 


04736 


;+03) 


-5 


11250 


12493.8 


1 


245 


+10) 


3 


745 


+10) 


4 


230 


;+io) 


8.731 


2 


04721 


;+03) 


-5 


11240 


12037.0 


1 


355 


+10) 


3 


819 


+10) 


4 


302 


;+io) 


8.728 


2 


04706 


;+03) 


-5 


11230 


11571.5 


1 


475 


+10) 


3 


908 


+10) 


4 


361 


;+io) 


8.725 


2 


04693 


;+03) 


-5 


11220 


11105.7 


1 


601 


+10) 


4 


025 


+10) 


4 


387 


;+io) 


8.723 


2 


04679 


'+03) 


-5 


11210 


10813.9 


1 


687 


+10) 


4 


111 


+10) 


4 


384 


;+io) 


8.720 


2 


04666 


'+03) 


-5 


11191 


10698.9 


1 


735 


+10) 


4 


141 


+10) 


4 


375 


;+io) 


8.719 


2 


04641 


;+03) 


-5 


11171 


10596.2 


1 


779 


+10) 


4 


167 


+10) 


4 


368 


'+10) 


8.718 


2 


04616 


;+03) 


-5 


11132 


10452.9 


1 


849 


+10) 


4 


199 


+10) 


4 


357 


'+10) 


8.716 


2 


04567 


;+03) 


-5 


11053 


10280.9 


1 


954 


+10) 


4 


228 


+10) 


4 


341 


'+io) 


8.710 


2 


04470 


^+03) 


-5 


10975 


10109.0 


2 


060 


+10) 


4 


254 


+10) 


4 


333 


;+io) 


8.705 


2 


04375 


;+03) 


-5 


10881 


9903.3 


2 


189 


+10) 


4 


282 


+10) 


4 


332 


;+io) 


8.698 


2 


04264 


'+03) 


-5 


10788 


9699.1 


2 


321 


+10) 


4 


306 


+10) 


4 


338 


'+10) 


8.691 


2 


04156 


;+03) 


-5 


10600 


9288.2 


2 


598 


+10) 


4 


349 


+10) 


4 


365 


;+io) 


8.678 


2 


03946 


;+03) 


-5 


10413 


9128.8 


2 


769 


+10) 


4 


348 


+10) 


4 


361 


;+io) 


8.661 


2 


03742 


'+03) 


-5 


10126 


8883.6 


3 


032 


+10) 


4 


350 


+10) 


4 


361 


;+io) 


8.637 


2 


03438 


'+03) 


-5 


09864 


8658.1 


3 


276 


+10) 


4 


356 


+10) 


4 


367 


;+io) 


8.614 


2 


03168 


'+03) 


-5 


09502 


8480.8 


3 


533 


+10) 


4 


349 


+10) 


4 


351 


;+io) 


8.593 


2 


02804 


;+03) 


-5 


08913 


8365.4 


3 


831 


+10) 


4 


320 


+10) 


4 


331 


;+io) 


8.571 


2 


02225 


'+03) 


-5 


08323 


8250.0 


4 


126 


+10) 


4 


298 


+10) 


4 


309 


'+10) 


8.550 


2 


01657 


;+03) 


-5 


07604 


8125.0 


4 


476 


+10) 


4 


274 


+10) 


4 


285 


;+io) 


8.525 


2 


00979 


;+03) 


-5 


07245 


8062.5 


4 


654 


+10) 


4 


263 


+10) 


4 


274 


;+io) 


8.512 


2 


00645 


'+03) 


-5 


06885 


8000.0 


4 


832 


+10) 


4 


254 


+10) 


4 


265 


;+io) 


8.500 


2 


00314 


;+03) 


-5 


05586 


7870.0 


5 


399 


+10) 


4 


217 


+10) 


4 


228 


;+io) 


8.460 


1 


99143 


;+03) 


-5 


04287 


7740.0 


5 


978 


+10) 


4 


188 


+10) 


4 


200 


'+10) 


8.420 


1 


98002 


;+03) 


-5 


02682 


7635.0 


6 


690 


+10) 


4 


173 


+10) 


4 


186 


'+10) 


8.320 


1 


96630 


;+03) 


-5 


01077 


7530.0 


7 


422 


+10) 


4 


173 


+10) 


4 


186 


;+io) 


8.220 


1 


95298 


;+03) 


-4 


98543 


7420.3 


8 


517 


+10) 


4 


192 


+10) 


4 


206 


;+io) 


8.099 


1 


93263 


;+03) 


-4 


96121 


7322.3 


9 


622 


+10) 


4 


233 


+10) 


4 


249 


;+io) 


7.976 


1 


91371 


'+03) 


-4 


93694 


7247.4 


1 


079 


+11) 


4 


295 


+10) 


4 


311 


;+io) 


7.834 


1 


89535 


;+03) 


-4 


88839 


7116.9 


1 


330 


+11) 


4 


452 


+10) 


4 


471 


;+io) 


7.577 


1 


86012 


;+03) 


-4 


82271 


7005.9 


1 


704 


+11) 


4 


687 


+10) 


4 


709 


;+io) 


7.323 


1 


81490 


;+03) 


-4 


76027 


6900.2 


2 


138 


+11) 


4 


918 


+10) 


4 


943 


;+io) 


7.081 


1 


77394 


;+03) 


-4 


66481 


6786.4 


2 


962 


+11) 


5 


261 


+10) 


5 


293 


;+io) 


6.733 


1 


71493 


'+03) 


-4 


55600 


6698.4 


4 


214 


+11) 


5 


595 


+10) 


5 


636 


;+io) 


6.357 


1 


65188 


;+03) 


-4 


43506 


6600.7 


6 


164 


+11) 


5 


857 


+10) 


5 


909 


;+io) 


5.938 


1 


58639 


;+03) 


-4 


30139 


6535.7 


9 


154 


+11) 


6 


080 


+10) 


6 


152 


'+10) 


5.523 


1 


51861 


;+03) 


-4 


15567 


6475.9 


1 


383 


+ 12) 


6 


411 


+10) 


6 


509 


'+10) 


5.096 


1 


44911 


;+03) 


-4 


07753 


6443.9 


1 


716 


+ 12) 


6 


619 


+10) 


6 


735 


'+10) 


4.877 


1 


41303 


'+03) 


-3 


99939 


6411.8 


2 


123 


+12) 


6 


888 


+10) 


7 


025 


;+io) 


4.659 


1 


37795 


;+03) 


-3 


91534 


6380.7 


2 


660 


+12) 


7 


303 


+10) 


7 


468 


;+io) 


4.432 


1 


34124 


;+03) 


-3 


83129 


6355.5 


3 


318 


+12) 


7 


839 


+10) 


8 


037 


'+10) 


4.219 


1 


30545 


;+03) 


-3 


65298 


6256.8 


5 


318 


+12) 


8 


810 


+10) 


9 


096 


;+io) 


3.777 


1 


23355 


'+03) 


-3 


47354 


6127.5 


8 


539 


+12) 


9 


151 


+10) 


9 


567 


;+io) 


3.356 


1 


16501 


'+03) 


-3 


30000 


6135.4 


1 


313 


+13) 


1 


141 


+11) 


1 


203 


;+n) 


2.960 


1 


10119 


;+03) 


-3 


13428 


6112.0 


1 


982 


+13) 


1 


376 


+11) 


1 


465 


:+n) 


2.592 


1 


04187 


;+03) 


-2 


97768 


6044.7 


2 


936 


+13) 


1 


514 


+11) 


1 


637 


;+n) 


2.262 


9 


87399 


;+02) 
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-2.83040 


5955 





4 
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+13) 


1 


566 


,+H) 


i 
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+11) 


1 


961 


9.37651(+02) 


-2.68992 


5851 
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6 


070 


+13) 


1 


550 


:+n) 


i 


759 


+11) 


1 


689 


8.91585(+02) 


-2.55407 


5736 


4 


8 


528 


+13) 


1 


467 


:+n) 


i 


722 


+11) 


1 


538 


8.48236(+02) 


-2.42406 


5589 


2 


1 


187 


+14) 


1 


274 


:+n) 


i 


569 


+11) 


1 


416 


8.07902(+02) 


-2.29839 


5444 


7 


1 


635 


+ 14) 


1 


076 


:+n) 


i 


410 


+11) 


1 


299 


7.70062(+02) 


-2.17917 


5291 


1 


2 


222 


+ 14) 


8 


586 


'+10) 


i 


230 


+11) 


1 


193 


7.35240(+02) 


-2.07437 


5132 


7 


2 


923 


+ 14) 


6 


353 


;+io) 


i 


047 


+11) 


1 


108 


7.05562(+02) 


-2.02936 


5061 


5 


3 


291 


+ 14) 


5 


088 


'+10) 


9 


455 


+10) 


1 


073 


6.93045(+02) 


-1.98435 


4987 


5 


3 


708 


+ 14) 


4 


206 


;+io) 


8 


838 


+10) 


1 


041 


6.80717(+02) 


-1.90597 


4858 


6 


4 


566 


+ 14) 


3 


061 


'+io) 


8 


291 


+10) 





985 


6.59768(+02) 


-1.83566 


4746 


6 


5 


502 


+ 14) 


2 


052 


;+io) 


8 


032 


+10) 





937 


6.41444(+02) 


-1.77067 


4646 


6 


6 


534 


+ 14) 


1 


359 


'+10) 


8 


206 


+10) 





894 


6.24891(+02) 


-1.70766 


4549 


5 


7 


723 


+ 14) 


8 


764 


;+09) 


8 


713 


+10) 





851 


6.09191(+02) 


-1.64196 


4469 


2 


9 


155 


+ 14) 


5 


754 


;+09) 


9 


576 


+10) 





809 


5.93154(+02) 


-1.56959 


4402 


8 


1 


099 


+ 15) 


4 


027 


;+09) 


1 


081 


+11) 





763 


5.75801(+02) 


-1.48889 


4328 


7 


1 


348 


+ 15) 


2 


794 


;+09) 


1 


235 


+11) 





713 


5.56793(+02) 


-1.39929 


4297 


7 


1 


670 


+ 15) 


2 


267 


'+09) 


1 


451 


+11) 





669 


5.35970(+02) 


-1.30098 


4280 


3 


2 


105 


+ 15) 


2 


136 


;+09) 


1 


742 


+11) 





626 


5.13289(+02) 


-1.18674 


4277 


6 


2 


742 


+15) 


2 


223 


'+09) 


2 


168 


+11) 





584 


4.87053(+02) 


-1.05256 


4295 


8 


3 


722 


+15) 


2 


624 


'+09) 


2 


837 


+11) 





540 


4.56269(+02) 


-0.90681 


4354 


9 


5 


137 


+15) 


3 


876 


'+09) 


3 


911 


+11) 





527 


4.22570(+02) 


-0.76297 


4423 


9 


7 


043 


+15) 


5 


966 


;+09) 


5 


406 


+11) 





524 


3.88814(+02) 


-0.62093 


4492 





9 


620 


+15) 


8 


962 


'+09) 


7 


427 


+11) 





520 


3.54964(+02) 


-0.48037 


4586 


4 


1 


302 


+16) 


1 


527 


'+10) 


1 


034 


+12) 





555 


3.20827(+02) 


-0.34126 


4683 


4 


1 


755 


+16) 


2 


637 


'+1°) 


1 


434 


+12) 





593 


2.86306(+02) 


-0.20325 


4782 


3 


2 


359 


+16) 


4 


675 


'+1°) 


1 


982 


+12) 





637 


2.51306(+02) 


-0.06709 
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